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SOME  TOPICS  ON  LIMITERS  AND  FM  DEMODULATORS 


by 


Maung  Gyi 


Systems  Theory  Laboratory 
Stanford  Electronics  Laboratories 
Stanford  University  Stanford,  California 


ABSTRACT 


An  analysis  is  marie  of  wide-deviation  frequency-modulation  systems 
having  low  nominal-carrier-frequency  to  information-bandwidth  ratios. 

Since  limiting  plays  an  important  role  in  such  systems,  the  effects  of 
hard  limiting  of  many  signals  in  random  noise  are  analyzed.  Expressions 
are  given  for  the  output  signal- to-signal  ratios  (SSR),  the  output  signal- 
to-noise  ratios  (SNR),  and  the  power  in  any  crossproduct.  The  effect  of 
the  power  in  each  output  component  as  a  function  of  the  input  SNR  is 
investigated.  It  is  found  that  for  all  values  of  input  SNR's  greater 
than  10,  the  strengths  of  the  various  output  components  are  relatively 
constant.  For  the  case  of  more  than  two  input  signals,  a  weak  signal¬ 
boosting  effect  manifests  itself  when  the  input  SNR's  are  less  than  1 
and  the  input  SSR's  are  greater  than  l/lO.  The  signal-suppression  effect 
and  the  signal- to-signal  power  sharing,  together  with  their  dependence  on 
input  signal  noise,  are  presented  for  various  cases. 

Expressions  are  presented  for  the  autocorrelation  function  and  the 
spectral  power  density  of  hard-limited  FM  pulse  trains,  which  allow  the 
computation  of  the  intermodulation  products  in  the  information  bandwidth 
or  in  any  other  frequency  interval.  The  effect  of  the  nominal-carrier 
frequency  on  the  interference  ratio  is  exhibited  graphically  for  various 
constant  deviations. 

A  partial  interaction  of  the  positive  and  negative  spectra  of  the 
modulated  wave  causes  extraneous  outputs.  For  all  such  spurious  components 
to  be  filterable,  a  relation  is  derived  between  the  carrier  frequency,  the 
maximum  frequency  of  information,  and  the  spectral  bandwidth  of  the  modu¬ 
lated  wave.  The  analysis  results  in  systems  which  have  better  performance 
capabilities . 
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PRINCIPAL  SYMBOLS 


k,m, m„. .m 

12  p 


E(x) 


f 

c 


f 

m 


(N/2)‘'/^  h  „  „  [see  Eq.  (2.29)] 

k  .  .nij  '• 

r5(cos  ilf(t)]  *  5(f  +  f^) 

3(cos  ilr(t)]  *  6(f  -  f^) 

i  [1  +  sgn  (f)] 
c 

^  [1  -  sgn  (f)] 
c 

complete  elliptic  integral  of  the  second  kind 

„  /  Ji\  r,  2  .  2  ,1/2  ^ 

=  E  ^x,  —j  =  J  [1  -  X  sin  41]  '  d(D 

nominal-carrier  frequency 
modulating  frequency 


f _  maximum  frequency  present  in  the  information 

mm 

3{G]  Fourier  treuisform  of  [o] 


g(t)  *  g(t) 

I  (x) 
m 

IM(p ,n) 
IR(p,n) 


k  self-convolutions  of  g(t) 

modified  Bessel  function  of  the  first  kind  of  order  m 
(imaginary  argument) 

intermodulation  product  at  frequency  (2pf^  +  *^^m^ 

interference  ratio  at  frequency  (2pf^  + 

Bessel  function  of  the  first  kind  of  order  m  (real 
argument ) 


k 


percent  deviation 
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R^fO.t] 


two-dimensional  joint  characteristic  function 

probability  distribution  function  of  ( y) 

E[x(t  )  •  x(t  )]  =  autocorrelation  function  of  the  real 
X  2 

random  process  x(t) 

R  (t) 

X 

spectral  power  density  =  Fourier  transform  of  R(t) 

R('i)  dr,  u'  =  2, if 

lW 


S^(f) 

S(f)  • 

■X  [1  +  sgn  (f)1 

Sjf) 

S(f)  • 

■1  [1  -  sgn  (f)] 

^f 

[sin 

\If(t))  x  6(f  -  f^) 

c 

x(t) 

Hilbert  transform  of  x(t 

^,(t) 


1  f'xis 

”  J-  ,.  ^  - 


A  Af  /  k  c 

modulation  index  =  +  —  =  1+  *  y~ 

m  \  m 


Dirac  delta  function  =  unit  impulse  at  t  =  0 

jl  if  m  =  0 

Neumann  numbers  =  \ 

b  if  m  >  0 

random  process  defined  by  ['l'(t)  -  \|f(t  -  r)] 


R  ( X ) 

j  =  normalized  autocorrelation  function 


f)  e  “  df ,  a-  =  2nf 
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fui 
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a(f ) 


Fourier  transform  of 


-ia^T 

e 


dT 


cD  (f) 

y 


'Kt) 

lir(t) 


sgn  (f) 


characteristic  function  which 


E[e^“^] 


=/: 


p(y)  dy. 


is  equal  to 

OD  =  2nf 


information  in  the  case  of  phase  modulation 


information  in  the  case  of  frequency  modulation 


convolution;  that  is, 


g(t)  *  f(t)  =  J  g{a)  f(t  -  a)  da 


-r 


f(a)  g(t  - 


for  f  <  0 
for  f  =  0 
for  f  >  0 
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I .  INTRODUCTION 


The  purpose  of  this  research  is  to  analyze  some  of  the  problems 
peculiar  to  wide-deviation  frequency-modulation  systems  having  low  nominal- 
carrier  frequency  to  information-bandwidth  ratios.  Some  form  of  frequency 
modulation  is  employed  in  many  tape-recorder  systems  due  to  the  fact  that 
the  effect  of  random  envelope  variations  occurring  during  playback  can  be 
eliminated  at  the  receiver  by  hard  limiting  prior  to  demodulation.  A 
typical  electromagnetic  recording  channel  has  a  bandwidth  from  an  upper 
limit  of  a  few  megacycles  down  to  a  lower  limit  of  a  few  hundred  cycles. 

If,  .'f'or  example,  it  is  desired  to  record  information  down  to  dc,  some  form 
of  modulation  must  be  used  before  the  information  can  be  so  recorded.  The 
parameters  of  a  typical  channel  are  a  low  nominal-carrier  frequency  and 
the  widest  possible  deviation  ratio.  The  high  deviation  is  needed  to 
obtain  an  acceptable  signal- to-noise  ratio  (SNR).  Since  the  channel  band¬ 
width  is  fixed,  most  of  the  modulation  schemes  result  in  reduction  of  the 
information  bandwidth  that  can  be  recorded. 

In  order  to  maximize  ^he  information  bandwidth,  a  vestigial  technique 
is  sometimes  employed.  At  other  times  conventional  frequency  modulation 
is  used.  In  either  case,  spurious  outputs  arise  because  the  bandwidth  of 
the  channel  is  not  wide  enough  to  record  the  entire  portion  of  the  frequency 
modulated  spectrum  that  carries  significant  amounts  of  energy.  The  ultimate 
performance  of  such  a  system  is  often  limited  by  these  spurious  outputs 
rather  than  by  the  SNR.  In  other  instances,  many  signals  are  multiplexed 
on  the  same  channel.  Then  it  is  of  importance  to  be  able  to  compute  the 
strengths  of  various  crossproducts  and  the  signal-to-noise  ratios  of  each 
signal . 

A .  BACKGROUND 

The  effects  of  hard  limiting  on  signals  and  noise  have  been  studied  by 
various  investigators  from  time  to  time.  Using  the  techniques  developed 
by  Bennett  and  Rice  [Ref.  l]  and  by  Middleton  [Refs.  2,3],  Davenport 
studied  the  case  of  a  single  sine  wave  and  additive  gaussian  noise  when 
passed  through  a  v  law'  device  and  gave  explicit  results  for  the  changes 
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in  signal- to-noise  ratio  produced  by  an  ideal  limiter  [Ref.  4].  Later 
Granlund  [Ref.  5]  and  then  Baghdady  [Ref.  6]  studied  the  problem  of  hard 
limiting  of  two  sine  waves  without  noise  and  their  results  included  some 
expressions  for  the  evaluation  of  the  strengths  of  output  signals  and 
intermodulation  products.  Davenport  and  Root  [Ref.  7]  studied  the  prob¬ 
lem  of  bandpass  limiting  of  two  sine  waves  in  the  presence  of  additive 
gaussian  noise.  Their  investigations  were  reinforced  recently  by  Jones 
[Ref.  8],  who  studied  the  problem  in  great  detail  in  order  to  examine  the 
effects  of  limiters  in  long  pulse  radars  and  other  cases  where  substantial 
pulse  overlaps  occurred.  Jones  studied  the  case  of  two  signals  in  noise 
and  gave  relations  among  the  output  signals,  noise,  and  intermodulation 
products  due  to  bandpass  hard  limiting  in  that  case  only. 

The  present  study  extends  the  analysis  to  the  case  of  bandpass  hard 
limiting  of  p  signals  in  additive  gaussian  noise.  It  is,  in  short,  a 
generalization  of  the  problem  investigated  by  Jones,  using  the  nonlinear 
transform  technique  which  is  convenient  and  mathematically  rigorous.^ 

The  results  are  shown  to  agree  with  those  of  (l)  Jones,  when  p  =  2; 

(2)  Granlund  and  Baghdady,  when  p  =  2  and  when  the  input  signal-to-noise 
ratio  is  made  arbitrarily  large  (no-noise  case);  and  (3)  Davenport,  when 

p  =  1. 

B.  OUTLINE  OF  PRESENTATION 

Chapter  II  presents  an  analysis  of  the  effect  of  hard  limiting  on 
P  signals  in  noise.  Expressions  are  given  for  the  spectral  power  density 
and  the  autocorrelation  of  the  limited  output,  the  output  SNR's,  the  out¬ 
put  signal- to-signal  ratios  (SSR),  the  signal-suppression  effect,  and  the 
intensity  of  any  crossproduct.  Some  curves  and  numerical  results  are  also 
presented,  and  some  asymptotic  values  are  derived. 

Chapter  III  presents  a  derivation  of  expressions  giving  the  output 
speocrum  of  wideband  frequency-modulated  pulse  trains  generated  from 

^Prior  to  publication  of  this  study,  P.  D.  Shaft,  who  has  also  investigated 
the  problem,  presented  a  paper  at  the  IEEE  Convention  in  New  York  on  24 
March  1965  [Ref.  9].  Mr.  Shaft  was  formerly  with  Western  Development 
Laboratories,  Philco  Corporation,  and  is  now  with  the  Stanford  Research 
Institute,  Menlo  Park,  California. 


hard- limi ted  frequency-modulated  signals.  Expressions  are  given  for  the 
spurious  outputs  in  the  demodulated  signal.  A  system  of  indexing  the 
spurious  components  is  developed,  and  curves  are  given  which  show  the 
dependence  of  spurious  outputs  on  carrier  frequency,  deviation,  and  the 
information  frequency.  Interference  ratios  are  defined,  *hen  computed 
and  tabulated  together  with  experimentally  measured  values. 

Chapter  IV  presents  an  analysis  of  the  "aliasing"  phenomenon  which  is 
due  to  the  interaction  of  the  positive  and  the  negative  spectra  when  the 
frequency  of  the  nominal  frequency-modulated  carrier  has  a  value  lower  then 
half  the  spectral  bandwidth  of  the  modulated  wave.  Expressions  are  given 
that  relate  the  nominal-carrier  frequency  to  the  maximum  frequency  of 
information.  Analysis  and  block  diagrams  of  proposed  systems  are  pre¬ 
sented  which  are  improvements  over  existing  systems. 
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II.  HARD  LIMITING  OF  p  SIGNALS  IN  RANDOM  NOISE 


A.  FORMULATION  OF  THE  PROBLEM 

The  effects  of  hard  limiting  of  p  signals  in  random  noise  are  ana¬ 
lyzed,  assuming  the  input  conditions  described  in  the  following  five  sub¬ 
sections.  The  system  studied  is  shown  in  the  block  diagram  of  Fig.  1,  in 
which  it  is  assumed  that  the  input  signals  and  noise  are  bandlimited  by 

bandpass  filtering  centered  around  f  . 

c 


FIG.  1.  BLOCK  DIAGRAM  OF  THE  RECEIVER  ANALYZED. 

1 .  Limiter  Input  Random  Process 

The  limiter  input  random  process,  defined  as  x(t),  consists  of 
p  unrelated  sine  waves  and  additive  gaussian  noise.  In  particular,  for 
every  i  in  the  (l,p)  interval, 

s^(t)  =  y/2^  COS  (cD^t  +  0^)  (2.1) 

where  the  are  statistically  independent  random  variables,  each  having 

a  uniform  distribution  on  the  interval  (0,2n). 

The  input  x(t)  is  therefore 

x(t)  =  s^(t)  +  s^Ct)  +  ...  +  Sp(t)  +  n(t) 

P 

=  ^  cos  (o'^t  +  (!)^)  +  n(t)  (2.2) 

i=l 
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2 .  Input  Noise 


The  input  noise  n(t)  is  a  sample  function  of  zero-mean,  station¬ 
ary,  narrowband  gaussian  noise.  It  is  assumed  symmetrical  around  f^ 
since  the  signals  and  noise  are  bandpass  filtered  prior  to  limiting. 

3 .  Autocorrelation  Function  of  the  Input  Noise 

The  noise  n(t)  is  a  sample  function  of  a  stationary  random 
process.  The  covariance  function  of  n(t)  can  be  readily  written  as 
[Davenport  and  Root,  Ref.  7,  p.  169] 

R(t)  =  R^(t)  cos  a^^t  -  R^g(t)  sin  u.^t 


=  [x^t)  * 


1/2 


COS  a.  t  +  tan 
c 


,  R  (t) 

■1  cs 

TJT) 


(2.3) 


where  a.  =  2nf  . 

c  c 


Since  f  is  chosen  such  that  S(f)  is  symmetric  about  it. 


and 


s  (f  +  f  )  =  s  (-f  +  f^; 

+'  c  +  c' 


(2.4) 


S_(f  -  f^)  =  S_(-f  -  f^; 


(2.5) 


which  gives 


S  (f)  =  2S  (f  +  f  ) 
c  +  c 


S  (f)  =  0 
cs 


This  equation  implies  that 


(2.6) 


r( t )  =  2R  ( t )  cos  a  t 
'  '  c  c 


(2.7) 
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Adopting  the  usual  normalization  by  defining  p(t)  such  that  it  contains 
unit  power,  one  can  write  the  covariance  function  of  the  input  noise  as 

R  (t)  =  Np(t)  cos  (JO  T  (2.8) 

n  c 

4 .  Limiter  Transfer  Function 

The  limiter  output  y(t)  can  be  written  in  terms  of  its  input  as 


r 

for 

x(t) 

>  0 

y(t)  =  g[x(t)]  =  <  0 

for 

x(t) 

=  0 

(2.9) 

for 

x(t) 

<  0 

5 .  Narrowband  Assumption 

The  narrowband  assumption  implies  that  the  difference  between  any 

pair  [cu.iCJJ.]  chosen  from  the  set  {oo- ,(jd_ , .  .  .  ,ao  ,U3  }  is  negligible 
1  J  1  A  p  c 

magnitudewise  when  compared  to  any  member  in  the  set;  that  is, 

|<x^  -  ujj  I  «  for  any  i,j  f  (k)  (2.I0) 

where  {k]  =  {l,2,...,p,c).  This  assumption  guarantees  that  the  inter¬ 
ference  between  adjacent  spectral  zones  at  the  output  of  the  limiter  will 
be  negligible.  Hence  the  first  spectral  zone  can  be  obtained  by  suitable 
bandpass  filtering  around  f^. 

B.  LIMITER  OUlPUT 

1 .  Derivation  of  Autocorrelation  Function 

The  signal- to-noise  ratio  among  the  output  components  of  the 
limiter  output  y(t)  can  be  obtained  from  the  autocorrelation  function 
of  the  limiter  output  which  is  defined  as 

R  (t^)  =E^{g[x(t)]  g[x(t  +  t)])  (2.11) 

y  ^ 
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or 


R  =  E  (g[x(t)]  g[x(t')]}  (2.12) 

y  ^ 

where  t'  ^  t  +  T.  Using  the  techniques  suggested  by  Bennett  and  Rice 
[Ref.  l],  Middleton  [Ref.  2],  and  Davenport  and  Root  [Ref.  7],  one  can 
represent  the  transfer  function  of  the  nonlinear  device  by 


g(x) 


(2.13) 


where 


C  =  the  contour  w  =  u'  +  jv,  u'  >  u 

+  ^ 

C  =  the  contour  w  =  u"  +  Jv,  u"  < 


Hence 


g(x) 


e 


xw  dw 
w 


The  contours  and  C  are  shown  in  Fig.  2. 

It  should  be  noted  that  both  contours  include  the  entire  w  =  jv 
axis  except  the  origin,  which  is  excluded  because  of  the  presence  of  a 
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.;\i'  -.(•/  v'.v.v.v.'v.v,:/', 


*  *  PtiNE. 


FIG.  2.  CONTOURS  OF  INTEGRATION. 

pole  at  that  point.  The  integrands  are  identical  and  therefore,  for 
compactness,  the  equation  of  g(x)  is  represented  by  a  single  integral 
over  the  contour  C.  For  final  evaluation,  one  can  simply  evaluate  over 
and  C  separately  and  add. 


g(x) 


1  I  wx  dw 
2nj  ®  w 


(2.14) 


Therefore , 


"  ■; - (2  I  ~  I  H [w  x(t)  +  W  x(t')]] 

(2.15) 

By  definition  E^{exp  [w^x(t)  +  w^,x(t')]}  is  the  two-dimensional  joint 

characteristic  function  and  can  be  denoted  by  M  (w  ,w  .).  Hence 

X  t  t ' 


Ry(t,t') 


1 

(23lj  )^ 


(2.16) 


Since  the  input  process  has  an  independent  nature,  manipulation  of  M 

X 

is  particularly  convenient.  The  input  x( t )  is  the  sum  of  the  zero-mean 
statistically  independent  random  processes , 
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(2.17) 


P 

x(t)  =  n(t)  +  ^  s^(t) 
i=l 

and  therefore 

P+1 

=  n  \ 

i=l  1 

P+1 

^x^  {exp  [w^x^(t)  +  w^,x^(t')]}  (2.18) 

It  can  be  shown  that  the  Joint  characteristic  function  for  stationary 
gaussian  noise  is  [Ref.  7,  Eq .  (13.43),  p.  289;  Parzen,  Ref.  10] 

[I  (\  ^  '"t’)  "  Vf] 

Each  of  the  sine-wave  signals 

s .  =  J2S.  cos  (cc .  t  +  <t> .) 

has  a  joint  characteristic  function 

(Wt.Wt,)  =  E  [exp  (wt/is"  cos  0^  +  w^,  y2S~  cos  6^1  )]  (2.20a) 

i 

The  exponentials  are  expanded  using  the  Jacobi-Anger  formula  [Watson, 

Ref.  11;  also  see  Eqs.  (B.lOa)  -  (B.lOd)  on  p.  115,  this  report]  to  give 

00  00 

M  (w.  ,w  ,  )  =  V  \  e  c  I  (w  ^28.  )  I  (w  ,  ^28.  )  E(cos  me  cos  ne  ,  ) 
s.  t  t'  Z_,  Zj  m  n  m  tv  i'  n  t'^  x'  t  t' 

^  m=0  n=0 

(2.20b) 
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ff  viv><  knt  y-n  v>  yn  v>.  Ln  VM  yxv>i>  y>  yn  ->  LXir*  v?* 


The  expectation  can  be  simplified  to 


E(cos  me  cos  nO  ,  ) 

L  X* 


E[cos  m(cjD^t  + 


cos  nfou.  t '  +  <i> .)] 
1  1 


[(m  +  n)ou.t  +  naj.T  + 

1  1  1 


+ 


+  E  cos  [(m  -  n)to^t  -  nco^T  +  (m  -  n)(t)^|| 


when  n  ^ 


m 


cos  (niD.T) 

M  i 


when  m  =  n 


(2.21) 


The  above  simplification  gives 


=  2!  V  <='■*  (2-22) 


m=0 


where 


Ijjj(x)  =  [modified  Bessel  function  of 
the  first  kind  of  order  m] 


f  =  Neumann  numbers 
m 


2 

V 


if  m  =  0 

if  m  >  0 


From  Eq.  (2.19) 


exp  [ ( T )  w 


t't'i  =  2 


i  [R^(c)  WjWj,]'' 


k=0 
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Therefore, 


=  exp  (I  wj)  exp  (|  exp  fR_^(T) 


=  exp  (f  J)  exp  (I  wj.)  ^  w^.^, 

k=0 


(2.23) 


which  gives 


f  dw  r  dw 
R  (t,t')  =  - -  I  -  I  -  M  (w  ,w  ,  ) 

>'  (2X1)2  Jc  ‘t  *t'  *  ‘  ' 


(2xj) 


i_  (•  %  r 

,.,2  Jc  »,  Jc  "f 


C2xj) 


P 

n  ”s 

i=l  i 


exp  (I  .2)  exp  (|  wj,)  ^  i  fR^(T) 


k=0 


Simplifying, 


-  ,  X  1  Y  r  k-1  /N  2\ 

y  "  "t)  ""t 


i=l 


I  V*  (1  "^)  “"'f 


(2.24) 
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Eiwitw¥«innin<yiruvuwwv>w’'<VN  i.'N  uxvxvh  »>  v*  v 


I 

'i 

\ 


Define; 


P  i=l  1 


dw  (2.25) 


With  this  definition,  Eq.  (2.24)  can  be  rewritten  as 


Ry(t,t')  =  Ry(T) 


°°  ,  , 

R  (t) 


=  I  I  ■  2  ^ 


k=0  m^=0 


1,  I  M  ^  cos  m.CD  T 

k,m, ..mlJ  m.  li 

1  P\i=l  1 


which  gives  the  autocorrelation  function  of  the  limiter  output  as 


00  00 


R  (t) 
y 


=  I  I  ■  2  ^ 


k=0  mj^=0 


\  m  »  I  n  m,  CD.  T 

k ,  m  .  .  m  \  1 1  m  i  i 

1  P  \i=l  i 


(2.26) 


A  product  of  p  cosines  can  be  expressed  as  a  sum  of  cosines 


T  r 

1.  [‘=°®  (“ifi 


^  “202 


...  +  m  f;  ) 
p  p 


+  cos  (m.t  +  „  +  ...-me.  )  +  cos  (m  P  + 

11  £  £  PP  11  2  ^ 


+  •  ■  .  “  m  E  ,  +  m  e 


p-1  p-1  p  p^ 


+  cos  (in,  ,  +  ^  ■  •  ■  “  -  m  0  )  .  .  .  2*^  terms] 

II  2  2  p-1  p-1  p  p'  J 


[cos  1  ...  1  y.^)] 


(2.27) 
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where  the  final  bracketed  term  represents  the  sum  of  2^  ^  cosines  by 
definition . 


a.  Simplification  of  the  Contour  Integral  h. 

k  ,  m,  .  .  m 
1  P 

In  evaluating  the  contour  integral  h.  ,  the  first 

k ,m, . .m 
1  P 

question  to  resolve  concerns  the  pole  at  the  origin. 


k  ,m,  .  .m 
1  P 


1  f  k-1 


n 


m 

i=i  i 


exp 


/N  2\ 

(2 " ) 


dw 


if  k-1 
—  I  w 
2nj 


P 

n 

i=l  1 


f]  I  (./S') 


exp 


(1*^) 


dw 


(2.28) 


Note  that  for  w  -►  0,  w^/2^p!  (see  Watson,  Ref.  11 ),  the  above 

integrand  becomes 


(k+Im.-l) 

w 

- — —  as  w  >  0 

constant 


The  integrand  is  well  behaved  along  the  whole  of  the  w  =  jv 
axis  except  possibly  at  the  origin.  The  pole  at  the  origin  will  also 
vanish  if  either  k  or  any  one  of  the  m^  is  greater  than  or  equal  to 
1.  If  the  dc  term  is  eliminated  from  the  covariance  function,  no  pole 
will  occur  at  the  origin.  All  the  limiter  inputs  have  zero  mean,  making 
the  integrand  analytic  along  the  entire  w  =  jv  axis.  The  two  contour 
integrals  can  now  be  reduced  to  a  single  integral.  Hence, 
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Converting  to  a  real  integral  by  using  the  relations  I 

r 

and  w  =  jv  yields 


-  —  [  /  N  2\  ^'"i  ij 


,  (k+Zm.-l)  -O)  p 

7  j  /  “ '  f  ^ )  n 

i=l  i 


^  (k+Zm^-1^ 


n 


j  v"-'  exp  (-  I  v^)  JI  J 

i  =1 


-  J  exp  (-  I  v^)  [j  J 

i=l 


J  (k+Im.-l) 


(k+£m.  ) 

1  -  (-1) 


Now 


JT  v"-*  exp  (-  I  v')  ^11 


(k+Zm. ) 

1  -  (-1) 


2  for  k  +  >  m.  odd 


= 


0  for 
V 


k  * 


even 
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(jz)  =  j*"  J^(z) 


(v  y 2S .  )  j  dv 


(v  y2S^  )  dv 


^  (vy2S.  )  dv 
i 

(v  y/2S .  )  dv 
i 


(v/is")  dv 


and  hence 


(k+Z.ni,-l) 


h  =  < 

k  .  .tn  \ 

1  P 


‘  J  v'^'^  exp  (-  I  v^)  ;I  (v/isT)  dv  for  /k  + 

■'0  1=1  I  V  / 

I-.) 


odd 


for  Ik  + 


(2.29) 


b. 


Computation  of  the  Integral  ^ 

’  1  ■  ■  p 

It  is  convenient  to  define  a  quantity 


(2.30) 


In  the  autocorrelation  function,  this  quantity  appears  with  an  index  of 
2,  hence  the  sign  of  the  term  will  be  dropped  and  the  quantity  redefined 
as 


fl 


k-1 


V 


exp 


k ,  m 


.  m 


P 


= 


(2.31) 


0 


otherwise 


It  has  been  found  preferable  to  adopt  the  following  two  forms; 

1.  Strong  Noise  Case.  N  is  greater  than  for  all  i  r  (l,p).  By 

an  appropriate  change  in  variable  and  subsequent  simplification, 


(2.32) 


15 
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It  is  easy  to  express  these  forms  in  terms  of  an  infinite 
series  involving  confluent,  gaussian,  or  generalized  hypergeometric  func¬ 
tions.  But  it  is  better  to  deal  directly  with  Eqs.  (2.3l)  -  (2.33)  when 
deriving  asymptotic  properties  and  limiting  values.  Due  care  must  be 
given  to  convergence  problems . 

For  numerical  results,  it  was  found  to  be  most  convenient  to 
carry  out  numerical  integration  on  the  computer  directly  rather  than 
resorting  to  the  evaluation  by  summation  of  hypergeometric  series."^ 


See  Appendix  A  and  Refs.  12-18  for  additional  formulas  and  tables  per¬ 
taining  to  integrals  involving  Bessel  functions. 
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c.  Wideband  and  Narrowband  Considerations 


Based  on  the  simplification  and  computational  procedures 


discussed  in  Secs,  la  and  lb  above,  the  following  expressions  that  are 


valid  for  any  bandwidth  are  obtained; 


00  00 


(T)  =  y  y  ■  ■  ■  y  . 

^  ^  Z-t  kl  k,m  .  .m  i  m. 

m  m  —A  ^  i  —  1  ^ 


COS  m.uj.T 
1  1 


k=0  m,  =0  m  =0 
1  P 


(2.35a) 


where  R  (t)  is  the  autocorrelation  function  of  the  input  noise. 


Using  the  expressions  and  notation  of  Eq.  (2.27),  one  can 


rewrite  the  above  equation  as 


uo  00 


P  ...  € 

m,  m 

1  P 


k=0  m^=0 


•"i  •  • ->P"  ^ 

Ip  2 


•  COS  [(m^^coj^ 


.  .  ±  m  OD  )t 

P  P 


(2.35b) 


The  Fourier  transform  of  Eq.  (2.35b)  gives  the  power  spectral 


density  of  the  limiter  output. 


00  e  .  .  .  e 

„  m,  m 


S  (f) 

y 


y  •••  y  — — [s  (f) s  ^f)] 

ZL  Z-  Z^  Vt  oP  k,m  .  .m  L  n'  '  n  '  J 

I.  K.4  J.p 


k=0  m^=0 


*  (S[f  -  (m^f^  ±  ...  ±  nipfp)]  +  6[f  +  (m^f^  ±  ...  ±  "1^1^)]) 

(2.36) 

where  |^S^(f)  *  implies  k  convolutions  of  with  Itself 


.  .  .  ±  m 


and  the  power  spectral  density  of  the  input  noise. 
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With  the  following  points  in  mind,  it  is  only  a  matter  of 
straightforward  computation  to  obtain  numerical  results  for  specific 
problems : 

1.  As  k  increases,  the  k  self-convolutions  of  S  (f)  may 

n 

rapidly  become  gaussian. 

2.  ^  reproducing  density,  e.g.,  exponential,  gaussian, 
etc . 

3.  Convolutions  with  delta  functions  are  particularly  easy. 

4.  In  practice,  one  is  usually  interested  in  a  certain  spectral  zone 
and  in  some  cases  in  a  few  terms,  which  simplifies  matters. 

5.  Use  should  be  made  of  the  fact  that  exists  only  for 

the  terms  where  k  +  Em^  is  odd  and  vanishes  for  the  rest.  Hence 
the  terms  to  be  computed  are  halved  at  the  outset. 

6.  The  value  of  ^k,mi  mp  small  as  to  be  insignificant. 

Up  to  this  point  no  use  has  been  made  of  the  narrowband 
restriction.  The  expression  for  the  autocorrelation  of  the  limiter  out¬ 
put  is  valid  for  any  bandwidth. 

Assuming  narrowband  restriction  and  normalizing, 

R^(t)  =  Np(t)  and 


where 


a(f)  =  r  p(t)  dT 


p(0)  J  o(f)  df  =  1 


SO  that  a(f)  is  a  probability  distribution  function.  The  input  noise 
has  a  density  which  is  even  about  f^, 

)  =  Np(-i )  cos 

SEL-65-056  -  18  - 


Substituting  in  Eq .  (2.35b) 


CX) 


=  I  I  -  S 


.  2 


c  .  .  .  e 
m,  m 


k=0  m^=0 


m  =0 
P 


1  P  ^ 


COS  [ 


±  n»  )t] 

P  P 


(2.37) 


To  simplify  the  above,  use  is  made  of  the  following  expression  [Ref.  7, 
page  298] : 


(k-l)/2 

S  (k-0!  1! 


05  T 

C 


for  k  odd 


i=0 


cos  (o5  T )  =/ 


k! 


k-1 


(k-2)/2 


I  oTTiTT: 


for  k  even 


(2.38) 


It  is  readily  seen  that  the  output  power  spectrum  is  centered  at  the  fol¬ 
lowing  carriers: 


(k  -  2i )  05  ±  mTa5.  ±  ...  ±  m  o) 

c  1  1  p  p 


(2.39) 
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Further 


I  I  I 


ks=l  m  =0 
k  odd 


bf  — i - E 

"'p=0  i  P  2  ^^0 


»'<”  I  irrrn-n 


cos  ([■(k  -  21  k  1  mo,  1  ...  +  mu.  It) 

c  11  P  p 


I  I  ■  I 


(  .  .  .  e 

m  m  , 

1  P  k 

-  P 


Ry(T)  =  < 


k=0  m  =0 
k  even 


k.m,  .m  D-1 

m  =0  1  P  2^ 

P 


(k-2)/2 

I  ~i]:  i: 


cos  {[(k  -  2i)cl  ±  ± 

c  11 


cos  [(ra^o,^  t  1  mpti.p)t] 


m  u  It}  +  - :: 

P  P  , «  2 


(10 


if  k 


if  k  + 


I-.) 


is  odd 


2 .  Spectral  Analysis  and  Filtering 

Since  the  autocorrelation  of  the  output  exists  for  only  the  odd 

integral  values  of  k  +  ,  the  spectrum  of  the  limiter  output  consists 

of  components  situated  in  zones  around  the  odd  harmonics  of  the  input 

frequencies.  As  stated  before,  the  narrowband  restriction  virtually 

eliminates  interference  of  these  output  zones  and  hence  any  one  zone  can 

be  selected  by  appropriately  bandpass  filtering  the  output. 

It  is  now  assumed  that  the  zone  around  the  first  harmonic  band 

around  f^  is  filtered.  The  filtered  output  process  is  denoted  by  z(t). 

The  general  expressions  for  R  {l)  and  S  (f)  are  cumbersome. 

z  z 

The  analysis  is  not  difficult  but  it  is  lengthy  and  requires  careful  manip 

ulations.  It  is  better  to  think  that  the  sums  of  cosines  constituting 

the  term  cos  (o )  are  the  carriers  which  are  shifted  up  or  down  by  the 

term  (11, -  .  .  .  ±  m  n  ),  To  get  the  f irst- harmonic  band,  the  following 
1  j.  p  p 

conditions  must  hold: 
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=  odd  integer  S  0 


(2.40) 


k  + 


(k  -  2i )  ±  m,  ±  . . .  ±  m  =  ±1 
1  P 


which  implies 


f 

0,2,4, .  .  . 

.  ,k+l  for 

k 

odd 

-H 

E 

.  .  .  +  m  1  =  < 
pi  ^ 

1,3,5, . . . 

,  ,k+l  for 

k 

even 

Now  for  each  of  the 

values  of  m.  which 

1 

satisfy  the 

above  equation,  a 

series  of  terms  is  written.  The  final  total  of  all  such  terms  in  the 
first- harmonic  band  may  in  some  cases  be  expressed  in  a  single  compact 
form.  This  is  possible  for  small  values  of  p,  such  as  0,  1,  or  2. 

For  instance,  this  procedure,  when  applied  to  the  case  p  =  2, 
gives  the  following  result  which  is  due  to  Jones  [Ref.  8,  Eq.  (13),  p. 
35] : 


■'zt-'*  ■  I 


I  s 


2  b 


’k|^|  M-|i+l|| 


k  +  i 


i=-M  k=  I  i  I  ,  I  i  I  +  2 


k  - 


P^(t) 


*  cos  [  i  o) 

c 


i  +  ijcWg  +  £{^2  -  )]t 


where 


k.m^m^ 


^(1) 


k.m^m^ 


Similar  expressions  can  be  derived  easily  for  p  =  0  or  1. 
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It  has  been  assumed  that  the  input  process  is  wide-sense  station¬ 
ary.  If  it  is  not  weakly  stationary,  the  only  change  that  would  result 
would  be  to  replace  T  by  (t  ,t  ). 

X  A 

3 .  Output  Signal- to-Noise  Terms 
a.  Output  Signal  Power 

Fortunately  it  is  not  necessary  to  calculate  the  series  for 
R  (t)  to  find  the  output  signal-to-noise  terms.  The  expression  for  the 

Z 

autocorrelation  of  the  output  consists  of  line  spectra  and  continuous 
noise  components.  The  periodic  part  of  the  output  is  primarily  the  result 
of  the  interaction  of  the  input  signals  with  themselves.  The  remaining 
terms  correspond  to  the  random  variations  of  the  output,  i.e.,  the  output 
noise.  These  terms  can  be  split  into  two  sets,  one  corresponding  to  the 
interaction  of  the  input  noise  with  itself  and  the  other  corresponding  to 
input  signals  interacting  with  the  input  noise.  Using  Davenport  and 
Root's  notation  [Ref.  7],  split  into  subsets: 


R,(T) 


=  R 


(t) 


I 

i=l 


R  (t; 

s .  n 
1 


P-1 

I 

i=l 


R  (t) 

s  .  s  ,  ,  n 
1  1  +  1 


R„„(t) 

nn 


(2.4i; 


where  Rs^..Sp(T) 
part  of  R^(t) ] , 
components . 


contains  the  output  signal  and  cross  terms  [periodic 
and  represents  the  direct  feedthrough  noise 


One  is  usually  interested  in  only  a  few  terms  belonging  to 
one  of  the  subsets.  Hence  it  is  not  necessary  to  calculate  the 

complete  sei’ies  of  R^(l),  which  is  usually  extremely  lengthy  for  values 
of  p  greater  than  1. 

1.  The  terms  in  Rg  ones  that  result  from  the  interaction 

i 

of  the  signal  i  with  the  input  noise. 

2.  The  terms  in  R„  „  n( r)  the  ones  that  result  from  the  inter- 

s  j  s  j  n 

action  of  both  tbe  input  signals  i  and  j  with  the  input  noise. 
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The  very  definition  of  these  terms  imposes  certain  obvious  conditions  on 
the  k  and  indices.  These  conditions,  together  with  the  need  that 
k  +  Em.  be  odd,  make  computing  the  terms  involved  a  simple  process. 

The  subscript  o  will  be  used  for  output  and  i  for  input. 
The  output  signal  power  is  given  by  imposing  the  conditions  k  =  0  and 
m.  =  0  for  j  ^  1,  which  give 


1  o 


th 
m . 

1 

zone 


2b 


2 

0,00. . m^O . . 0 


(2.42) 


Specifically,  for  the  first  zone,  m^  =  1  gives 


~  ^%,00.  .100.  .0 

first 

zone 


J  v"^  exp  I  v^j  J^Cv/is”) 


b.  Signal  and  Crossproduct  Components,  R  ,  ( ) 

S  .  .  b 

1  p 

One  can  easily  derive  expressions  for  signal  and  crossproduct 
componerts  in  both  the  wideband  and  narrowband  cases.  Inspection  of  the 
expressions  for  indicates  that  k  must  be  0.  Hence, 

in  the  wideband  unfiltered  case. 


R 

y  s 


1  ■ 


Y 

m^=0 


V 

m  :=0 
P 


^2 

0 ,m, . . m 
1  P 


m 

P 


cos 


.  .  .  i  m  (1 
P  P 


)t] 
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when  is  odd.  For  the  narrowband  case  this  reduces  to 


(  .  .  .  r 

m,  m 


R 

z  s 


(t)  =  y  ...  y  i — 

1--S  Z  Z  O.m  ..m  p-1 

^  m  =0  m  =0  ^  P  ^ 


cos  [  ±  .  .  .  +  m  ai  )  ] 

11  P  p 


with  the  additional  condition  imposed  on  the  indices  that 


m,  ±  .  .  .  1  m  =  ±1 
1  P 


For  p  =  1  the  above  reduces  to 


R  (t) 
y  s'  ' 


1 


mj^=l 


=  2 


s 


mj^=l 


b  cos  m,a)  'T 
0 ,  11 


=odd 


m^=odd 


which  is  precisely  the  result  given  by  Davenport  and  Root  [Ref.  7].  For 
the  first  zone,  obviously, 


zRsZ)  =  2bQ^^  cos  a.j 


For  any  bandwidth,  for  p  =  2, 


y  s, 

^12 


'Z  Z  ‘'0,m  m.,  "m/  m,  "’Zl  '  ">^2' 

m^=0  ^  ^  ^  ^ 

m  +m  =  odd 

X  4M 


For  the  narrowband  case  (first-harmonic  band),  ^R^  ^  (t),  the  additional 
condition  imposed  on  the  indices  is 


Z  S^Sg 
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m  ±  m  =  ±1 

X  ^ 


Hence , 


Sjf?, 


m  =0  m  =0 

X  A 

Zm.  =  odd 
1 


^'"2  "“l  "’2; 


COS  [  (m^u.^  +  m2(jj2)T] 


Further  simplification  in  this  case  is  possible  by  noting  that 
cos  [(m  oj,  +  m  aj_)T]  contributes  only  when 

XX  Ct  ^ 


mi  +  m^ 


=  1,  that  is, 


0,  -  1 

and 

0,  =  1 


/  2  2  \ 

This  gives  ^b^  cos  oj^^t  +  b^  cos  aigtj.  Also,  =  ±1 

implies  m  =  m  +  1.  This  contribution  to  the  first  zone  due  to  the 

M  X 

terms  arising  from  cos  [  (m  a:>  -  m  a.  )t]  is,  with  appropriate  restric- 

XX  A  ^ 

tions , 


1  II)  (fVi  - 

m,  =0  '  ^  ^  ^ 


=  y 


■"i  ^ 


^  y 


m^=0 


V  2^2 

,0 


“  ^^l  l"'rl|  ^  ‘’0,10  ‘"1^  '^0,01  “^2- 
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Hence, 


=  2  S  cos  {[m^Co^g  -  O)^)  -  o^glT] 

m^=-oo 


12 


(2.44) 

This  agrees  with  Jones'  result  [Ref.  8,  Eq.  (l5),  p.  36].  For  higher 
values  of  p,  the  computations  though  simple  and  straightforward  become 
lengthy.  Fortunately,  it  is  not  necessary  to  derive  the  above  series  for 
the  signal- to-noise  computations.  Additional  conditions  are  imposed  on 
the  indices  and  these  simplify  matters  extensively. 

c.  Direct  Feedthrough  Noise  Components,  R  (x) 

nn 

These  terms  impose  the  condition  that  m^  =  0  for  every  i. 
Since  the  autocorrelation  function  exists  only  when  the  sum  of  all  the 
Indices  is  odd,  k  must  be  odd.  Hence, 


R  (t) 

y  nn'  ' 


=  2 


\,00.  .0  k!  “c^ 


k=l 

odd 


00  (k-l)/2  ^2 


=  2  2 

k=l  i=0 

odd 


k  , 00 .  .  0  kl 

k : 


P  (k-i)l  i!  -  20^] 


(k-l)/2  2 


=  2  2 

k=l  i=0 

odd 


b -  p‘"(t) 


k  , 00 .  . 0  r  ,,  „  .  \ 

-  2i)V- 


For  (k  -  2i)  :=  1;  that  is,  i  =  (k  -  l)/2.  This  gives 


R  (x) 

z  nn 


k/  V 

P  (t) 


„  V  k,00.  .0  ^  ^ 

^  Z  kX  1 .  k  - 1  , 


k=l 

odd 
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The  expression  reduces  to  Eq,  (l7)  of  Jones  [Ref.  8]  when  p  =  2,  With 
no  restriction  placed  on  the  bandwidth,  these  terms  for  wideband  cases 
can  be  expressed  as 


R 

nn 


(t) 


I 

k=l 

odd 


r‘^(t) 

n 


k ,  00 .  .  0 


00 


h 


2 

k,00.  .0 


odd 


Also , 


S 

nn 


(f) 


k=l 


k,00.  .0 


odd 


C.  OUTPUT  FACTORS  DUE  TO  SIGNAL  AND  NOISE  INTERACTION 

These  computations  are  usually  lengthy  though  straightforward.  For 
example,  for  R„  n(’^)>  conditions  on  the  indices  which  must  be  met 

are: 

(i)  k  +  ^m^  =  odd  nonnegative  integer 

(ii)  m^  =0,  where  j  e  (l,p)  and  j  ^  i 

(iii)  m^  ^0  and  k  ^  0 
Hence  for  wideband  cases , 

OO  OO 

I 

k=l  m. =1 
1 

k+m.  =  odd 
1 
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R  ( 

s .  n' 
1 


2 

— — ; —  h, „  cos  m.cu.T 

k!  k,00. .m^OO. .0  i  i 


For  the  narrow  bandpass  case 


R  (r)  =  2 

y  s^n' 

I.V  1.V' 

I  s 

—  h^ 

k  !  k , 00 .  . m . 00 , . 0 

k  /  k 

p  (r)  cos  (.0  T  cos  m.Lij.T 

C  11 

k=l  m  =1 

1 

1 

k+m.  =  odd 

1 

(k-2)/2,  k  even 

tXi  iX' 

(k-l)/2,  k  odd 

=  2 

I  S 

k2  k. 

b  1  1  ' 

k,00.  .m.O.  .0  ' 

^  S  (k  -  j):  j! 

k=:l  m.=l 

1 

X 

j-o 

k+m,  =odd 

1 

• 

cos  { [ (k  - 

2j  )a^  -  m.u'.  ]  r  J 

Note  that  the 

cos  ( [ (k  - 

2 j  )a)  +  m.co.  ]'rj 

C  11^-^ 

terms  are  not  permissible 

owing  to  the  condition  on  the  indices. 


For  the  first- harmonic  band,  further  conditions  are  imposed, 


2j  -  m^  =  ±1 ; 


k  -  2j 


(m.  4  1) 


J  = 


k  -  (m.  ±  1) 
1 


k  -  j  = 


k  +  (m,  4  1 ) 
1  ' 


Also  when  k  is  even,  m^  =  odd,  which  implies  k  =  m.  ^  1.  Applying 
this  and  simplifying, 


R 

z  s 


=  s 


2b 


k,00. .  |i  +  l|0. .0 


k/  ^ 

p  (t) 


i=-cxi  k=  i  ,  i  +2 


k  + 


k 


cos 


(|i|a.^  -  |(i  +  l)|aj[ 
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This  result  reduces  to  that  of  Jones  [Ref.  8]  when  p  =  2.  Similarly, 
expressions  can  be  developed  for  the  other  factors  of  the  autocorrelation 
function  of  the  filtered  output. 


D.  SMOOTH  LIMITERS  WITH  ERROR- FUNCTION  TRANSFER  CHARACTERISTICS 

The  above  results  are  applicable  to  the  class  of  smooth  limiters 
which  have  a  transfer  function  describable  by  an  error  curve 


r 

g^(x)  =  erf  (x/a) 


g(x)  =  <0 


for 


for 


X  >  0 

X  =  0 


(2.45) 


^g^(-x)  for  X  <  0 


2v 

by  simply  replacing  N  in  the  above  analysis  by  (N  +  a  } ,  which  is 
seen  to  be  valid  for  any  value  of  p.  The  value  of  a  controls  the 
smoothness  of  the  limiter.  For  the  case  of  one  signal  in  noise,  see 
Galejs  [Ref.  19],  Blachman  [Ref.  20],  and  Jones  [Ref.  8], 


E .  ASYMPTOTIC  RESULTS 

1 .  Case  I,  Weak  Signals 


Since  all 
than  S^  for  all 


p  signals  are  buried  deep  in  noise, 
i  e  ( 1 , p ) .  From  Eq .  (2.3) 


N  is  greater 


^^0,00.  .  (mj^=l)0.  .0 


exp  (-X  /4) 
X 


dx 


(2.46) 
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For  further  evaluation,  use  is  made  of  the  following  arguments; 


1.  The  power  series  representing  the  Bessel  functions 


J 


(z) 


n=0 


is  absolutely  convergent  for  all  values  of  x,  real  or  complex, 
less  than  infinity. 

2.  Hence  the  power  series  obtained  by  the  multiplication  of  a  finite 
number  of  such  series  is  also  absolutely  convergent  for  all  finite 
values  of  x. 

For  detailed  proof  refer  to  Watson  [Ref.  11,  Sec.  5.41,  pp.  147-148].  In 
particular. 


where  ^2n  ’  ^  ^  ^  terms  involving  S^/n  and  its  higher  powers. 
This  series  is  absolutely  convergent.  Therefore  the  integration  and  sum¬ 
mation  can  be  interchanged,  giving 
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Since  each  integrand  exists  for  all  n  ^  0, 


Furthermore,  c^,  ,  etc.  involve  S^/N  and  its  higher  powers,  and 

since  N  can  be  arbitrarily  chosen  large  enough  such  that  the  higher 
order  terms  can  be  neglected,  then,  for  very  large  N  such  that  S^/N  -*•  0 
for  all  i  €  (l  ,p) , 


b 


0,00. .10, .0 


(2.47) 


Also,  under  these  conditions,  for  every  i  €  (l,p),  the  output  signal 
power  in  signal  i  in  the  first  zone  is 


(S.) 

1  o 


(2.48) 


Hence  under  strong  noise  conditions,  input  signal- to-signal  ratios  are 
prest.ved,  i.e., 

for  J,k  e  (l,p)  (2.49) 
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To  compute  the  signal- to-noise  ratio  at  the  output,  the  total 
output  noise  power  in  this  case  is  given  by  the  n  x  n  terms.  The  total 
noise  power  in  z(t),  defined  as  N^,  is 


N 

o 


=  R 

nn 


R  (0) 
s.n'  ' 
1 


+ 


s 


R  (0) 

s.s.  ,n  ' 
1  1+1 


+  R 


(0) 


From  Eq.  (2.32)  and  arguments  similar  to  the  ones  used  above, 


[  s/s./N  x] 


higher 

order 

terms 


dx 


As  N  -  00  such  that  S^/n  -  0,  then  ^  for  any  m.  g  1  -  0 

as  S^/N  -+  0.  Hence  b^^  ^  is  the  only  term  to  survive.  That  is, 

the  noise  at  the  output  is  contributed  by  n  \  n  terms  only  under  strong 
noise  conditions,  in  which  case, 


which  by  definition  of  the  gamma  function  becomes 


\  ,  00 .  .  0 


(2.50) 
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The  output  noise  is  given  by 


R  (O) 
z  nn 


=  2 


V  X  ^^(1) 

^  k  +  1  ,  k  -  1  , 


k=l 

odd 


Substituting  k  =  2J  +  1, 


00  ^2 


N  =  R  (O) 
o  z  nn 


A  V  ^  ^  2) 

„2  z  j:(J  +  1): 


j=0 


Using  the  generalized  factorial  notation  given  in  Ref.  13, 


R  (0) 
z  nn'  ' 


J=o 


ij  J 


-'2^1  (i  V  ■) 


=  —  for  N  -*  00  and  S .  /N  0 

2  i' 


(2.51) 


In  the  above  equation,  E(x)  denotes  the  complete  elliptic  integral  of 

the  second  kind.  This  result  could  have  been  predicted  by  observing  that 

the  limiter  output  is  a  (-l,+l)  square  wave.  Hence  the  total  power  in 

2 

the  fundamental  zone  is  (s/n  ).  Since  S.  -»  0  for  every  i  e  (l,p), 

2 

-»  (s/n  ).  When  N  -♦  uo  in  such  a  way  that  (S^/n)^  0,  the  output 

signal-to-noise  ratio  is  given  by 
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for  any  i  e  (l,p) 


(2.52) 


that  is,  a  loss  of  1  db  approximately.  Hence,  when  noise  controls  the 
limiter,  the  following  general  results  apply: 

1.  All  input  signal- to-signal  ratios  are  preserved  when  hard  bandpass 
limiting  occurs,  i.e.,  signal  suppression  effect  is  absent. 

2.  The  signal-to-noise  ratio  at  the  output  for  every  signal  is  1  db 
less  than  the  signal-to-noise  ratio  at  the  input. 

Particular  cases  of  the  above  have  been  derived  previously  by  Davenport 
for  the  one-signal- in-noise  case  [Ref.  4]  and  by  Jones  in  the  two- 
signals-in-noise  case  [Ref.  8]. 

2 .  Case  II,  One  Strong  Signal 

a.  Signal  and  Crossproduct  Power 

In  this  case,  one  signal  is  very  strong  relative  to  the  noise 
and  to  each  of  the  other  signals.  Let  the  strongest  signal  be  This 

is  no  restriction  since  the  numbering  of  the  input  signals  is  arbitrary. 
Here 


P 

(s )  s  y  s  +  N 

J.  1 

S.--2 

It  is  also  given  that 


/  \  2 

The  limiting  value  of  the  output  strong  signal  power  [S, j  is  2b  , 

X.  O  vy  I  ^  v/  •  «  vJ 

evaluated  under  the  above  given  conditions.  From  Eq .  (2.33), 


exp 


(-Ny^/4sJ 


Jj(y) 


i=:2 


(iy 


The  integrand  exists  for  all  values  of  positive  indices  on  y  and,  by 
the  arguments  developed  previously. 
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(S./S^)  -  0,  i  ^  1, 


and 


JmCO 

L  i  J^(y)  dy  for 


(S^/N) 


Since  the  integrand  is  1,  ^  2/jt  and  therefore, 


(^)o  -  (I)  = :! 


which  is,  as  expected,  the  total  energy  in  the  first-harmonic  zone. 
An  alternate  way  to  arrive  at  the  above  is 


Noting  that  ~  e^/y  2nz  for  large  z. 


'^0,10.  .0 


2 

n 


which  is  the  above  result. 

The  limiting  value  of  any  of  the  weak  signals  at  the  limiter 

output  is 


(S.  )  = 

1  o 


2b 


0,0. . (m^=l)0.  .0 


^  1 
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where 


j/i 


Since  ^  “*  easily  be  justified  by 

use  of  the  above  arguments  that 


0,0. 


.  (m.=l)0. .0 


JqCy)  dy 


This  integral  is  the  well-known  Sonine  and  Schafheitlin  integral.  Hence, 


0,0. . (m.=l)0. .0  n 


r  (I) 

2r(2)r  (i)  2  1 


i  2-  M 
2’  2’ 


ip 

n  \  S. 


1/2 


as  S ,  /S  -»  0 

i'  1 


(2.54) 


Therefore,  for  (S  /n).  oo  and  (S./S  )  0, 

1^1  1'  1  i 


(2.55) 
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Hence,  the  output  weak-signal  to  strong-signal  power  ratio  is 


for  (Sj/n).  -00  and  °  (2.56) 


Therefore,  if  the  strong  signal  controls  the  limiter,  then  all  the  other 

signals  are  modified  such  that  the  output  weak-signal  to  strong-signal 

power  ratio  is  one-fourth  of  the  input  power  ratio.  In  the  case  of  two 

signals  without  noise,  a  similar  result  is  given  [Ref.  8]. 

It  is  easy  to  evaluate  any  interfering  component  under  the 

above  conditions.  It  would  be  convenient  to  adopt  the  notation 

S  ^  ^  to  mean  the  crossproduct  component  at  a  frequency 

00.  .m^^O.  .m.O.  .0 

(m  0)  -  m  oj  ) .  An  example  would  be  the  evaluation  of  the  intermodulation 

i  i  J  j 

component  which  occurs  at  a  frequency  of  -  uj^ ) .  It  is  best  to 

start  with  the  expression  for  s  (r)  in  such  cases.  Therefore, 

for  (S^/N)^  -00  and  (S^/S^).  -0  for  i  •=■  (2,p),  the  ratio  of  the 
output  power  of  the  strongest  intermodulation  component  to  the  output 
power  of  the  weak  signal  is 


200. .iOO. .0 
S. 


-  1 


(2.57) 


When  p  signals  are  present  at  the  input  of  the  limiter, 
such  that  the  strongest  signal  controls  the  limiter  and  all  the 

others  are  weak  signals,  the  output  power  of  the  (2aij^  -  )  intermodu¬ 

lation  component  approaches  that  of  the  power  of  the  weaker  signals 
[Ref.  8]. 

b.  Noise  Power 

To  determine  which  of  the  terms  in  the  output  power  contribute 
to  the  noise,  the  following  analysis  is  adopted: 


N 

o 


R 

z  nn 


(0) 


+ 


R 

z  s j^n 


(0) 


+ 


R 

z  s^n 


(0) 


+ 


R 

z  s. 


s  n 
P 


(0) 


-  37  - 


SEL-65-056 


Using  Eq.  (2.33),  it  can  easily  be  shown  that 


P 

1  I 

]  ! 

i=2 


m./2  m. 
(s^/s^)  "  y  ' 


m. 

2  ^  m,  ! 
1 


1 


higher 

order  (y)  dy 

terms  1 


Since  S 


i/^1 


0  if 


m. 


i  /  1, 


b,  0 

k , m,  .  .  m 

1  p 


if  any  /  0,  i  t,  (2,p) 


(2,58) 


The  above  equation  implies  that  the  terms  which  contribute  to  the  output 
noise  power  are  n  ■  n  terms  and  s^  x  n  terms  and  therefore 


k.m^OO. 


J  (y)  dy 


Since  (S^/n)^  ->  oo  and  ^  terms  involving  the  lowest 

power  in  k,  i.e.,  k  =  1  are  retained.  The  first-order  component 
belonging  to  the  n  '  n  subset  is 


2b 


2 

1,00- -0  „ 

i:  o: 


2 

2 

n 


(2.59) 


The  rest  of  the  terms  in  the  n  n  subset  containing  higher  powers  of 

(n/Sj^)^  are  neglected.  Also,  irrespective  of  the  value  of  m^ ,  which 

must  be  even  for  b,  „  to  exist, 

1  .mj^OO.  .0  ’ 
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1 

n 


1/2 


1 ,m^00. .0 


"  fe)  i  S*’'’ 


Furthermore , 


OO  OO 


R  „(o) 


y 


k=l 


2  k 

h,  nr,  n  COS  U  T  COS  m  03,  T 
k.m^OO. .0  c  11 


T=0 


R  „(o) 

z  s^n 


=  2  (I  h^,200..o)  2  "OS  03^'^  "OS  2a3^T 


k=l, 

m^=2 


T=0 

filtered 


=  2b 


1 , 200 . . 0 


2_  /JN 
2  \S 


(2.60) 


Hence  for  (S^^/n)^  -*  oo  and  ^  ^  (2,p),  the  total  noise 

output  power  is 


2  2  4  /  N 

^o  ""  ^*^1,00.  .0  ^*’1,200.  .0  ~2  is 

n  \  1/ 


The  output  signal-to-noise  ratio  for  the  strong  signal  is 


(2.61) 


The  3-db  increase  in  the  output  SNR,  which  occurs  due  to  hard 
bandpass  limiting,  applies  only  to  the  strongest  signal,  which  controls 
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the  limiter.  This  result  for  the  case  of  one  signal  was  given  by 
Davenport  [Ref.  4],  and  this  analysis  also  checks  with  the  one  by  Jones 
[Ref.  8]. 

For  wealt  signals  ,  i  ^  1, 


(2.62) 


Each  of  the  remaining  (p-l)  weak  signals  suffers  an  SNR 
loss  of  3  db  when  the  limiter  is  controlled  by  the  strong  signal. 


3 .  Case  III,  p  Signals  of  Equal  Strength 


a.  Signal  Power 


From  Eq.  (2.43)  it 
depends  solely  on  bi.  „  _  . 

K  ,  .  .  lUp 

equal  strength  gives 


is  seen  that  the  output  signal  power 
The  condition  that  all  the  signals  are  of 


^0 , 00 . . ( 


m.=l)0. .0 


Ji(.v)  Jq  ^(y)  dy 

(2.63) 


As  a  result  of  the  symmetry  of  the  integrand,  it  can  be  concluded  that 
all  signal  output  powers  are  equal; 


o 


exp 


2 


J’i(y)  Jq  ^(y)  dy 


For  the  noiseless  case  and  p  =1  or  2,  closed-form  solutions 
exist  [Refs.  7,8].  For  all  other  cases,  solutions  can  be  developed  in 
the  form  of  an  infinite  series  involving  generalized  hypergeometric  func¬ 
tions  . 

For  large  p,  the  following  procedure  culminates  in  a  useful 
approximation . 
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Define 


ine  l(p)  =  j  y  ^  exp  Jj(y)  Jq  ^(y)  dy 


Comparing  the  series 


with 


2  4 

Jo(y)  =  1  -  ^  +  V— 

2  2*4 


/  2  4 

j  (v)  -  ^  _ +  y _ 

1'"’^'  2\  2*4  2*4*4*6 


exp 


-(I)  --^^(1)  IT 


exp 


=  1- 


2  •  4 


and  u^^suming  p  large,  one  obtains  in  the  region  of  interest, 


Hence , 
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By  a  suitable  change  of  variable, 


Since  the  integration  has  been  carried  out  to  infinity,  it  is  preferable 
to  replace  [p  -  l/2  +  (n/S^^)]  in  the  denominator  by  [p  +  (n/S^)], 
giving 


The  normalized  output  signal  power  is  by  definition 


(8/„2) 


1.05  +  10 


in  decibels  (2.65) 


This  equation  implies  that  the  output  signal  power  in  decibels  decreases 
linearly  with  the  number  of  signals.  Also,  if  the  noise  power  is  equal 
to  the  signal  power,  then  the  output  signal  power  for  p  signals  with 
noise  is  approximately  the  same  as  for  (p+l)  signals  without  noise. 

To  verify  the  above  results  and  also  to  justify  the  validity 
of  the  approximation,  the  exact  output  signal  power  was  computed  by 
numerical  integration  on  a  Philco  2000  computer  and  compared  with  the 
results  obtained  from  the  above  approximate  formula.  The  above  formula 
gives  results  which  are  within  0.1  db  of  the  exact  value  for  p  greater 
than  or  equal  to  4.  The  results  are  tabulated  for  values  of  p  as  high 
as  100  (see  Table  3,  p.  61 ) . 


b.  Crossproduct  Power 


The  output  power  of, 

interference  frequency  of  |m  a) 

'  i  i 


^0. .m^O. .mj .  .0- 

-  m  .ou  .  I  ,  is  given  by 
J  J 


which  has  an 
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appropriate 

term 


«  «  =  R  (o) 

0. .m.O. .m^ . .0  z  s,s„..s 
1  J  1  2  p 


=  2  r  y-*  J 

L"  Jo  ”i 


l2 


(y)  (y)  Jq  ^(y)  dy 

j 


=  s 


0. .m .0, .m. . .0 
J  1 


(2.66) 


which  by  definition  has  an  interfering  frequency  of  |m.a^^  -  |  .  Also 


note  that 


m  .cu . 

J  J 


m  .  =m .  - 1 
J  1 


|m.a)i 


m  .o) .  +  CO  .  I 

J  J  J ' 


00.  +  CO  , 

_i _ J 

2 


+ 


and 


I  Vi 


m.co , 

1  J 


m  ^  =m . - 1 

J  1 


m.cjo.  -  CO.  -  m.a' , 

11  1  1  J 


00.  +  CO  . 

_ J 

2 


It  can  be  concluded,  therefore,  that  the  output  crossproduct  terms  situated 
on  either  side  of  the  arithmetic  mean  frequency  have  the  same  power. 
Similarly,  it  can  be  shown  that  the  enitre  output  spectrum  in  the  first 
zone  is  symmetrical  about  the  arithmetic  mean  frequency. 

Other  results  can  also  be  derived  when  the  inputs  have  equal 
powers ,  e . g . , 
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"power 

at 

1  m .  a' . 

-  m  .ai 

1  1  1 

J  q 

power 

at 

1  m  .a' . 

-  m.a  . 

Output  power  at 

Im.o).  -  m.oj.l  =  ( 

'  J  1 

1  J 

J  J '  \ 

1 

power 

at 

m  ,ui . 

-  m.a' 

'  J  1 

1  q 

etc . 

Therefore,  if  some  of  the  signals  are  of  equal  strength,  the  number  of 
components  to  be  computed  decreases  very  rapidly. 

For  the  case  p  =  2,  it  is  easy  to  compute  the  entire  series 
comprising  the  signal  outputs  in  the  first  zone.  The  formula  (see  Ref. 
21) 
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which  is,  as  it  should  be,  the  total  energy  in  the  first  zone.  This 
special  result,  which  states  that  the  power  contained  in  the  intermodula¬ 
tion  components  falls  off  as  the  reciprocals  of  the  squares  of  odd  inte¬ 
gers,  was  given  by  Jones  [Ref.  8].  The  output  power  in  the  crossproduct 
=11.. 100.. 0  is  ei™  by 


q  p-q 
times  times 


iQoP-q 


norm 


iV^ 

_8_ 

2 


(2.68) 


The  integrand  is  best  evaluated  by  direct  integration  on  a  computer.  To 
find  an  approximation  to  the  above  integral  for  large  values  of  p,  the 
following  procedure  was  adopted  Equations  (2.64)  give 


With  an  appropriate  change  of  variable  and  subsequent  simplification 


Il(p) 


'o  /  q 


2  dz 


p-i*s 


N 


-  z 
e 
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Since  the  integration  has  been  carried  out  to  infinity,  it  is 
ferable  to  replace  (p  -  q/2  +  N/S^^)  in  the  denominator  by 
1/2  +  N/Sj^)  to  give 


I^(p) 


which  gives 


S 

iQoP-q 


lOq  logj^Q 


in  decibels 


Equation  (2.69)  gives  the  following  results  in  decibels: 
For  q  =  1 : 


^®100. .0^0 

^  — 

r  /  N^ 

1.05  +  10  log  fp  .  - 

1 

norm 

L  \  ly 

For  q  =  3: 


/  \1 

^®11100. .0^0 

- 

7.07  .  30  log  /p  -  1  +  Aj 

norm 

L  \  i/J 

For  q  =  5: 

r  /  \n 

^®111110. .0^0 

a:  - 

3.55  +  50  log  (p  -  2  +  |- ) 

norm 

L  \  l/J 

This  result  agrees  with  Eq. 


(2.65)  . 
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again  pre- 
p  -  q/2  + 


(2.69a) 


^  (2.69b) 


(2.70a)^ 


(2.70b) 


(2.70c) 


( 2 . 70d  ) 


For 


For 


q 


«  “ 

/  N^ 

90  log  ( p  -  4  +  - 

1  -  15.29 

norm 

L  \  1/ 

- 

(2.70e) 


All  the  above  expressions  give  results  which  are  very  close  to  those 

obtained  by  numerical  integration  on  the  Philco  computer  for  large  values 

of  p.  Some  of  the  results  are  given  in  Tables  5,  6,  and  7,  pages  63-65. 

The  output  power  in  the  crossproduct  ^  „  a  r  given  by 

2  1^0^"^ 


2riqoP-q-r 


norm 


f 


-1  ^ 

y  dy  exp 


>;(y: 


p-q-r 
0 


(y) 


A 


(2.71) 


Proceeding  on  similar  lines  and  approximating  the  Bessel  functions  in  the 
region  of  interest  as  indicated  by  Eqs .  (2.64)  results  in 


By  a  suitable  change  of  variable  and  simplification, 


I2(P) 


.r+1 


P  - 


q 

2 


(2r+q)/2 


(^) 
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Hence , 


(2.72a) 


(2.72b) 


If  r  =  0,  then  Eq .  (2.72b)  reduces  to  Eq .  (2.69) 
If  r=l,  q=l,  then 


norm 


13.09  +  30  log^Q  / p  -  I  +  |- 


in  decibels  (2.73) 


If  r  =  1,  q  =  3,  then 
S 


norm 


9.57  +  50  log^^  -  I  ,  A 


in  decibels  (2.74) 


If  r  =  1,  q  =  5,  then 
S 


2HV-® 


norm 


1.61.  70  p-l-f 


in  decibels  (2.75) 


If  r=2,  q=l,  then 
S 


2  1  p-3 
2  10*^ 


norm 


15.59  .  50  log^  /p  -  I  .  f 


in  decibels  (2.76, 


The  above  expressions  give  values  which  are  close  to  those  obtained  by 
integrating  numerically  on  the  computer.  Some  of  the  results  are  shown 
in  Tables  5  and  7  on  pages  63  and  65. 
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It  may  be  added  that  such  techniques  can  be  used  to  get 
approximate  expressions  for  the  case  when  the  signal  strengths  are  not 
equal . 


F .  COMPUTED  RESULTS 

For  p  ?'  2,  Jones  used  the  technique  of  expanding  the  integrand  in 
a  series  involving  hypergeometric  functions,  integrating,  and  using  the 
resulting  series  for  computational  purposes  [Ref.  8].  This  procedure 
was  found  to  be  unsuitable  for  values  of  p  >  2.  All  computations  were 
performed  on  the  Philco  2000  computer  by  carrying  out  the  appropriate 
integrations  numerically.  The  numerical  integration  technique  was  also 
used  by  Shaft  [Ref.  9].  Some  of  the  results  presented  here  follow  the 
same  format  as  those  in  Ref.  9. 

The  following  results  are  presented  in  tabular  and  graphical  form;'!' 


Figures  3a  and  3b  show  the  output  power  in  the  various  signal 
components  for  various  p  =  3  and 

when 

a.  One  input  signal  is  varied  such  that  S^  '  (S^  =  S^). 

b.  Two  input  signal.s  are  varied  such  that  =  S^)  '  S^. 

Figures  4a-4c  show  the  effect  on  the  output  power  in  the  various 
signal  components  when  p  =  4  and  when 

a.  One  input  signal  is  varied  keeping  the  other  (p  -  l)  =3 
signals  equal  and  fixed. 


Two  input  signals, 
S 


and 


'2’ 


are  varied,  keeping  S. 


and 


constant  in  such  a  way  that  (S  =  S 
4  1  ^ 

Three  input  signals  are  varied  in  such  a  manner  that 

-  Sj 


(S3  =  S^) 


5. 


^®1  "2  “3'  ~4' 

Figures  5a- 5c  show  the  power  in  the  output  signal  components  for 

the  case  of  p  =  5,  under  the  conditions  specified  therein. 

Tables  1  and  2  give  the  signal-suppression  effect  in  the  case  of 

p  =  3  and  p  =  4 . 

Figures  6a-6b  and  7a-7c  show  the  signal-suppression  effect  graphi¬ 
cally  . 


^Figures  3-9  and  Tables  1-7  appear  at  the  end  of  this  chapter  beginning 
on  p .  54 . 
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6.  Table  3  gives  the  output  signal  power  in  the  case  of  p  equal 
signals  (no  noise)  for  values  of  1  :  p  100.  Both  computed  and 
approximate  values  using  Eq .  (2.65)  are  given. 

7.  Table  4  gives  the  output  signal  power  in  the  case  of  p  equal 
signals  in  noise  for  3  p  :  100. 

8.  Table  5  gives  the  strength  of  some  of  the  crossproducts  for  various 
input  signal-to-noise  ratios  in  the  case  of  four  equal  signals. 

9.  Table  6  gives  the  crossproducts  for  p  equal  signals  in  noise  for 
large  p. 

10.  Table  7  gives  other  crossproducts  for  the  case  of  p  equal  signals 
(no  noise) . 

11.  Figures  8  and  9  display  the  output  power  in  the  crossproducts  for 
large  p  in  the  equal-signal  case  and  clearly  show  that  the  values 
obtained  by  the  approximate  formulas  are  indiscernible  from  the 
exact  values  obtained  through  computer  evaluation  of  the  exact 
expressions . 

G.  CONCLUSIONS  TO  CHAPTER  II 

It  has  been  shown  that  many  of  the  limiting  properties  for  the  case 
of  hard  limiting  of  p  inputs  in  random  noise  are  the  same  as  those  for 
the  case  of  p  =  2  under  similar  conditions.  The  following  conclusions 
are  in  order; 

1.  The  limits  on  signal-to-noise  ratio  are: 

For  one  signal  in  noise, 


For  p  signals  in  noise, 
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2.  For  the  strong-noise  case,  the  following  results  which  are  known 
for  p  =  1,  2  [Ref.  8],  have  been  shown  to  be  applicable  for 
any  p ; 

a.  The  output  signal-to-signal  power  ratios  are  the  same  as  the 
input  signal-to-signal  power  ratios  for  any  number  of  inputs, 
that  is,  there  is  an  absence  of  signal  suppression. 

b.  For  any  p  the  signal- to-noise  ratio  at  the  output  is  1  db 
less  than  the  signal-to-noise  ratio  at  the  input. 

3.  For  one  strong  signal  and  (p-l)  weak  signals  S^: 

a.  If  the  input  power  ratio  is  very  large,  the  behavior 

approaches  that  of  the  one-signal- in-noise  case  as  the  ratio 
of  the  input  strong  signal  to  the  weak  signal  increases,  and 

P 

Si  »  ^  S.  .  N 

i=2 

Then,  for  large  input  signal-to-noise  ratios: 

(1)  The  range  of  the  normalized  SNR's  approaches  the  one- 
signal-in-noise  case. 

(2)  The  suppression  effect  manifests  itself,  with  the  amount 
approaching  the  limiting  value, 

(S,/S.)  -*4(S,/S.).,  that  is,  6  db. 

1'  1  o  '  1'  1 '1 

(3)  The  strong  signal  receives  an  improvement  in  relative 
SNR,  which  approaches  a  value  of  2,  that  is,  3  db. 

(4)  Each  of  the  (p-l)  signals  suffers  a  loss  in  SNR  by  a 
factor  of  2,  that  is,  3  db. 

b.  If  the  ratio  larger  than  approximately  20/l  and 

for  (Sj^/n)^  greater  than  1,  the  effects  of  limiting  are  the 

same  as  the  equivalent  two-signals-in-noise  case. 

c.  If  the  ratio  (S^/S^)^  is  greater  than  or  equal  to  20/l ,  and 

that  for  (S^/n).  is  less  than  1,  then  for  p  >  2,  the  sup¬ 

pression  effect  becomes  less  pronounced  and  ultimately  the  weak 
signals  are  boosted  Instead  of  being  suppressed. 
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4. 


5. 


6. 


7  , 

8. 


For  two  strong  signals,  say  S  and  S  ,  and  (p-2)  weak  signals 
S^,  the  behavior  follows  the  patterns  set  by  the  two-signals- in¬ 
noise  case.  There  is  one  exception,  however,  which  is  the  boosting 
of  the  weak  signals  for  (S  /S  ).  >  20/l  and  (S  /n).  <  1.  These 
are  the  very  same  conditions  as  for  item  3d.  The  strong  signals 
cause  a  beat  pattern  and  every  time  the  envelope  of  the  two  strong 
signals  falls  below  the  noise  or  weak  signal,  the  latter  gets  con¬ 
trol  of  the  limiter.  Consequently,  all  the  weak  signals  get  trans¬ 
mitted  without  any  attenuation.  The  effect  is  similar  to  the  one 
described  by  Jones  [Ref.  8],  which  is  an  attempt  to  give  some 
physical  explanation  as  to  the  decrease  in  the  relative  signal-to- 
noise  ratios  for  high  input  signal-to-noise  ratios.^ 

In  all  cases  the  crossproducts  follow  the  curves  in  Fig.  4  of  Jones 
[Ref.  8]  very  closely.  The  signal  power  in  the  crossproducts 
remains  nearly  constant  for  input  SNR's  greater  than  6  and  decreases 
very  rapidly  for  smaller  input  SNR's. 

From  the  numerical  results  and  graphs  presented,  it  is  concluded 
that  the  power  in  the  output  components  remains  relatively  inde¬ 
pendent  of  the  input  SNR's  for  input  strong-signal- to-noise  ratios 
greater  than  10  db. 

For  large  p,  and  approximately  equal  signals,  Eqs .  (2.70)  through 
(2.76)  are  valid. 

For  p  large,  the  following  results  were  obtained: 

a.  The  power  in  any  signal  or  crossproduct  term  decreases  in  a 

linear  fashion  when  plotted  on  a  log-log  basis  [see  Eqs.  (2.70)  - 
(2.76)]  . 


& 

\ 


\ 


f 

t 

i 


b. 


The  sum  of  the  powers  in  all  the  p  signals  at  the  output 
independent  of  p.  From  Eq.  (2.70), 


is 


1 

F, 

h: 


«  -1 .05  db  (2.77)  5 

V 

-  t 

"^Other  cases  may  be  reduced  in  a  similar  manner  to  predict  the  effects 
of  limiting.  ^ 

j 
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that  is,  the  power  contained  in  the  p  signals  is  approximately 
1.1  db  less  than  the  total  power  in  the  output  first- harmonic 
zone. 

c.  The  p-signals- in-noise  case  has  approximately  the  same  behavior 
due  to  hard  limiting  as  the  case  of  [p  +  (n/S^)]  signals  with¬ 
out  noise. 

d.  The  power  in  decibels  in  any  term  at  the  output  decreases  as 


-10 


log 


10 


(2.78a) 


that  is, 


Output  power  ~  p 


(2.78b) 
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FIG.  3.  OUTPUT  SIGNAL  POWER  FOR  THE  THREE-SIGNAL  CASE. 


TABLE  1.  SUPPRESSION  EFFECT  FOR  THE  THREE- SIGNAL  CASE  (p  =  3) 

(All  values  in  decibels) 
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TABLE  3.  OUTPUT  SIGNAL  POWER  FOR  p  EQUAL  SIGNALS  AND  NO  NOISE 

,  2 

Total  energy  in  first  zone  =  8/jt  =  0  db 
(All  computed  values*  in  minus  decibels) 


Number  of 
Signals 
(p) 

Exact 

Value 

Approximate 

Value 

Number  of 
Signals 

(p) 

Exact 

Value 

Approximate 

Value 

1 

0 

1.05 

14 

12.5 

12 . 5 

2 

3.9 

4.1 

15 

12.8 

12.8 

3 

5.6 

5.8 

16 

13.1 

13.1 

4 

6.9 

7.1 

17 

13.3 

13.4 

5 

7.9 

8.0 

18 

13.6 

13.6 

6 

8.7 

8.8 

19 

13.8 

13.8 

7 

9.4 

9.5 

20 

14.1 

14.1 

8 

10.1 

10.1 

30 

15.9 

15.8 

9 

10.6 

10.6 

40 

17.1 

17.1 

10 

11.0 

11.1 

50 

18.1 

18.0 

11 

11.5 

11.5 

75 

19.9 

19.8 

12 

11.8 

11.8 

100 

21.2 

21 . 1 

13 

12.2 

12.2 

* 

Values:  Exact 


-  (Sj 


1  O 


2 


Ji(y) 


Approximate  =  -[1.05  +  10  (p  +  N/S^)1, 

N  =  0,  noiseless  case. 

Note:  For  all  values  of  p  greater  than  4,  the  approximate  formula 

gives  the  answers  to  within  0.1  db.  For  the  most  part  the 
error  is  due  to  rounding  off  the  values  to  the  first  place. 
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TABLE  4.  OUTPUT  SIGNAL  POWER  FOR  p  EQUAL  SIGNALS  IN  NOISE 
(.All  computed  values*  in  minus  decibels) 


TABLE  5. 


CROSSPRODUCTS  VS  (s/n)  FOR  THE  FOUR- SIGNAL  CASE 


S,  =  S2  =  S3  =  S^ 


(All  computed  values  in  minus  decibels) 


Normalization  = 

00 

to 

II 

0 

db 

Notation;  =  S^^^3 

=  signal  at  frequency  If^^ 

^.100  -  =  O'— 

at  frequency  (2f^ 

'  '2) 

®2U1  =  %lj3 

=  power 

at  frequency  (2f^ 

.  .  f 

f  ) 

3  4' 

Crossproduct 

(s/N)i 

00 

20  db 

15  db 

10  db 

5  db 

0  db 

'iV 

6.9 

6.9 

6.9 

7.0 

7,2 

7,9 

28.5 

28.5 

28.6 

28.9 

29.7 

31.9 

21.5 

21.5 

21.6 

21.8 

22.5 

24.7 

Vl' 

31.6 

31.7 

32.1 

33.0 

35.9 

43.7 
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TABLE  6.  CROSSPRODUCTS  FOR  LARGE  p  IN  NOISE  (EQUAL-SIGNAL  CASE) 
(All  computed  values  in  minus  decibels) 


TABLE  7.  CROSSPRODUCTS  FOR  LARGE  p,  NO  NOISE  (EQUAL- SIGNALS  CASE) 


The  approximate  formulas  are  justified  for  large 


-(+«  ond+IOdb)-(S/N), 


Ok* 


20  30  40  50  60  70  80  90  100 

NUMBER  OF  SIGNALS,  p 


FIG.  8.  CROSSPRODUCT  S  „  „  VS  p  FOR 

iV"^ 


EQUAL- SIGNALS- IN-NOISE  CASE. 

In  Figs .  8  and  9  where  the  approximate  curve 
does  not  coincide  with  the  exact  curve,  it 
is  indicated  by  a  dashed  line. 
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Ill .  OUTPUT  SPECTRA  FOR  WIDEBAND  FM  PULSE  TRAINS 
GENERATED  FROM  HARD-LIMITED  FM  SIGNALS 


A.  INTRODUCTION 

In  most  narrowband  FM  systems  the  nominal  carrier  frequency  is  much 
larger  than  the  maximum  information  bandwidths,  e.g.,  an  FM  broadcast 
channel  bandwidth  is  approximately  36  kc,  while  the  carrier  is  a  few 
megacycles.  In  the  process  of  demodulation,  the  received  signal  is  sub¬ 
jected  to  hard  limiting.  The  amount  of  energy  in  the  sideband  of  the 

third  harmonic  which  lies  in  the  first- harmonic  band  is  typically  less 
-6 

than  10  of  the  energy  in  the  fundamental  band.  Hence  the  output  of 
the  limiter,  to  all  Intents  and  purposes,  consists  of  discrete  spectral 
zones  as  shown  in  the  following  sketch;  which  is  a  typical  spectrum  of  a 
limited  narrowband  FM  signal.  When  such  a  signal  is  passed  through  a 
narrow  bandpass  filter,  the  fundamental  band  can  be  filtered  out.  Since 
it  is  essential  that  such  a  bandpass  filter  have  an  excellent  time-delay 
characteristic,  the  filtered  fundamental  zone  does  not  suffer  any  appre¬ 
ciable  phase  distortion,  nor  is  the  spectral  distribution  of  the  first 
zone  disturbed.  As  a  result,  the  filtered  output  has  no  envelope  struc¬ 
ture.  Therefore  this  signal  can  be  demodulated  by  means  of  a  relatively 
narrowband  discriminator,  which  may  be  a  tuned  circuit  or  a  balanced 
detector  of  the  conventional  type  [Ref.  22].  Such  systems  have  been 
extensively  analyzed  by  various  investigators  [e.g.,  Middleton,  Ref.  2, 
and  Loughlin,  Ref.  23]. 
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As  an  illustration,  consider  a  case  of  monotone  frequency  modulation 
where  the  nomnalized  carrier  frequency  is  1  and  the  modulating  frequency 
is  approximately  0.5.  Assume  further  that  the  deviation  is  ±50  percent. 
The  limited  spectral  output  is  as  shown  in  Fig.  10a.  It  is  clear  that 
in  order  to  filter  out  the  fundamental  zone  immediately  prior  to  detec¬ 
tion,  as  before,  it  would  be  preferable  in  this  case  to  use  a  lowpass 
filter  of  bandwidth  2(Af  +  f^)  =  2  [Refs.  24  and  25]  .  The  filtered  out¬ 
put  will  then  be  as  shovm  in  Fig.  10b.  Because  of  the  high  selectivity 
of  the  lowpass  filter,  excessive  phase  distortions  in  the  sidebands  of 
the  filtered  output  will  occur.  In  addition,  since  the  filtered  output 
contains  quite  a  few  significant  sidebands  of  the  third-harmonic  zone, 
an  undesirable  envelope  structure  will  be  present.  The  detected  output 
may  be  computed  using  the  techniques  given  in  Refs.  24  and  26.  In  any 
case,  since  the  amplitudes  of  the  in-band  third-harmonic  lower  sidebands 
are  appreciable,  the  detected  output  will  certainly  have  unacceptable 
spurious  components.  Hence,  filtering  the  limited  output  is  not  desir¬ 
able  in  low-carrier  wide-dev iation  systems. 

Since  the  spurious  outputs  usually  limit  the  system  performance,  it 
has  been  necessary  to  define  an  "interference  ratio"  (see  Sec.  E).  A 
complete  analysis  of  such  a  demodulation  scheme  is  covered  in  this  chapter. 
The  autocorrelation  function  and  the  power  spectral  density  of  the  output 
process,  as  well  as  closed-form  expressions  for  the  interference  ratio, 
are  derived  for  the  general  case.  The  closed-form  expressions  are  used 
to  compute  the  spurious  outputs  which  are  within  2  db  of  the  experimen¬ 
tally  measured  values.  The  specifications  of  some  of  the  demodulators 
designed  on  the  basis  of  this  analysis  are  also  presented. 

B.  INPUT  PROCESS 
1 .  Assumption 

It  is  assumed  that  the  limiter  has  a  sufficiently  wide  bandwidth, 
which  would  imply  that  the  output  of  the  limiter  can  be  considered  to  be 

"^For  these  systems,  satisfactory  performance  may  be  obtained  through 
detection  by  means  of  integrating  a  train  of  arbitrary  but  fixed,  uni¬ 
directional  pulses  which  are  generated  every  time  the  limited  wave 
passes  through  zero  [Cady,  Ref.  27]. 
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switching  between  the  normalized  values  of  -1  to  +1  with  negligible  rise 
and  fall  times.  This  assumption  is  usually  valid  in  all  practical  systems. 


2 .  FM  Pulse  Train 

If  any  FM  signal  is  an  input  to  the  limiter,  then  the  train  of 

pulses  generated  by  placing  an  arbitrary  but  fixed  pulse  at  every  cross- 

t  ll 

over  or  at  every  p  crossover  (where  p  is  any  positive  integer)  of 
the  limited  signal  is  termed  an  "FM  pulse  train." 

An  arbitrary  but  fixed  pulse  implies  that  the  pulse  shape,  dura¬ 
tion,  amplitude,  etc.,  may  be  arbitrarily  chosen,  but  that  once  chosen, 
they  remain  invariant.  The  duration  of  the  pulse  is  such  that  the  suc¬ 
cessive  pulses  do  not  overlap  at  any  instant. 

It  may  be  noted  that  if  f^  is  the  nominal-carrier  frequency 

and  ±f  is  the  peak  frequency  deviation,  then  the  FM  train  with  a  pulse 
t  h 

at  every  p  crossover  has  a  pulse  repetition  rate  of  2f^/p  and  a 
pulse  deviation  of  ±2f/p. 

Conceptually,  and  in  practice  as  well,  one  of  the  best  ways  to 
generate  or  analyze  an  FM  train  is  indicated  in  the  block  diagram  of  Fig. 
11.  The  system  analyzed  is  described  in  the  following  sections. 


FIG.  11.  BLOCK  DIAGRAM  OF  THE  SYSTEM  ANALYZED. 
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a.  Input  Random  Process,  x(t) 

Since  the  results  for  phase  modulation  can  be  written  down 
directly  from  the  results  of  frequency  modulation,  it  will  suffice  to 
analyze  the  FM  case  only.  Hence  the  input  i-andom  process  is  assumed  to 
be  any  signal  of  the  form 

x(t)  =  cos  [oL'^t  +  +  '|f(t)]  (3.1) 

where  x(t)  is  a  process  consisting  of  a  sinusoidal  waveform  (of  constant 
amplitude)  whose  phase  is  the  sum  of  the  random  process  'lf(t)  and  the 
random  variable  If  'l'(t)  is  the  original  information,  then  this 

process  represents  phase  modulation.  In  the  FM  case,  the  derivative 
process  \lf(t)  represents  the  original  information,  that  is, 

'Kt)  =  f  s(0  (3.2) 


where  g(t)  is  the  original  information. 
Define 


0,  =a)t  +  (t)  +  k 

1  c  c  f 


(3.3) 


and 


0 


2 


T)  + 


(3.4) 


and  therefore 


x(t)  =  cos 
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c.  Pulse  Train,  z(t) 

The  FM  pulse  train  z(t)  consists  of  a  train  of  positive- 
going  (or  negative-going)  rectangular  pulses  of  constant  height  and  of 
constant  width  T  placed  at  every  crossover  of  the  limiter  output  process 
y(t).  As  shown  in  Fig.  11,  the  output  of  the  limiter  is  delayed  by  a 
fixed  amount,  say  T.  The  amount  of  delay  T  must  be  small  compared 
with  a  period  of  the  unmodulated  carrier  frequency.  This  delayed  output 
is  subtracted  from  y(t)  and  the  result  passed  through  a  zero  memory 
squarer.  This  gives  the  pulse  train  z(t),  where 

z(t)  =  [y(t)  -  y(t  -  t)]^ 

=  y^(t)  +  y^(t  -  T)  -  2y(t)  y(t  -  T) 
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Using  Eq .  (3.6)  gives 


’.(t)  =  2[1  -  y(t)  y(t  -  T)] 


(3.7) 


Hence  the  expression  for  [y(t)  y(t  -  T)]  will  give  the  expansion  of  the 
pulse  train  z(t),  the  modifications  being  of  a  trivial  nature. 

C.  EVALUATION  OF  [y(t)  y(t  -  T)] 

Equation  (3.7)  clearly  shows  the  necessity  of  evaluating 
[y('t)  y(t  -  T)].  With  Eqs.  (3.1)  -  (3.5),  it  can  be  shown  that 


y(t)y(t-T)=(i)  ^  2 


(-1) 


(n'l+m2~2)/2 


cos  m  cos  ni  £ 


m^=l  m^^l 
odd  odd 


"’l"'2 


(3.8a) 


~  (m  +m  -2)/2 

=  y  ">009) 

\  n  /  2  ^  11  22' 

X  M 


m  =1  m  =1 
X  A 

odd  odd 


+  cos  -  m^Gg)]  (3.8b) 


A 


-  ^  *^2  ^  ^4  ^  ^  ^2p  ^ 


(3.8c) 


t  hi 

where  C  represents  the  2p  spectral  zone.  From  Eq .  (3.8c)  it  is 
2p 

clear  that  the  odd  spectral  zones  do  not  exist. 


1.  Expression  for  the  0^*^  Spectral  Zone,  C, 


0 


The  0^^  spectral  zone  imposes  the  condition  on  the  indices  that 


mi  ^  >"2  =  ® 
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with  m  and  m  each  being  positive  odd  integers.  Hence  the  only 

1  2 

contribution  to  this  zone  from  the  term  cos  (m^G^^  -  ">2^2^ 


m, -1 


C  =  (i)  4  y  ^  -  cos  m  (e,  - 

0  \n/  2  Z-  „2  1  '  1 


mj^=l 

odd 


2 


cos 


mj^=l 

odd 


m. 


(3.9) 


From  Bromwich  [Ref.  28,  p.  371], 


2 


cos  m^^X 


lim 

a-»0 


cos  m^^X 

~2  2 
1  m  -  a 
m^^sl  1 


2 


n  sin 


=  lim 


(I-) 


na 


a-*0  4a  cos 


=  3(1-*) 


0  <  X  <  ji 


(3.10) 


Hence 


^0  "  ^  ■  f  ^®1  ■  ^2^  °  ^®1  ■  ®2^  ^ 

2.  Expression  for  the  2p^^  Spectral  Zone,  p  /  0 

The  contributions  to  the  2p^^  spectral  zone  can  be  written  down 
as  the  total  of  the  cos  (m^e^  +  m^e^)  terms  and  the  cos 
terms . 
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a.  The  cos  Terms 

These  terms  can  be  obtained  by  imposing  the  condition  on  the 

indices  that 


"*1  ^  "*2  ~  ™2  ”  "  "*1 

Since  m  and  m  can  each  take  on  only  positive  odd  integer  values 
1  2 

and  since  p  can  take  all  positive  integers  greater  than  zero,  only 
those  values  of  m^^  are  acceptable  which  satisfy  the  condition 

1  S  m^  S  (2p  1) 


The  application  of  this  condition  gives  the  subtotal  as 


=  (if  i  ^  'Vl  * 


(m^+2p-mj^-2)/2 


'2p 


m, 


mj^=l  ,odd 

j=(2p-m^) 


2  2p-l 

(I)  i  1 


m^=l 

odd 


m^(m^  -  2p) 


cos 


Og)  +  2pe2] 


(3.12) 


b.  The  cos  Terms 

These  terms  can  be  obtained  by  using  the  condition 


=  ±2p 


Consider  that  =  2p,  that  is,  =  m^^  -  2p .  Also  by  noting 

that  m  and  m  can  have  all  odd  integer  values  greater  than  zero  and 
1  2 

the  p  can  only  have  positive  nonzero  integer  values,  it  can  be  shown 
that 
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00 


(2p  +  1)  S  5 


Therefore 


2p 


(-1) 


[m^+(m^-2p)-2]/2 


mj^=2p+l  ,odd 


- J - — -r -  COS  [m  0  -  (m  - 


2p)e2] 


m. 


=  (in^-2p) 


-{€\  1 


mj^=2p+l 

all  odd 


f_l)P 

- 7-^ - '  ~—r  cos  [m,  (e, 

m^(m^  -  2p)  1^  1 


O2)  +  2pe2] 


(3.13) 


Now  consider  (m^  -  m2)  =  -2p,  which  implies  m2  =  (m^  +  2p).  Imposing 
the  following  conditions, 


p  g  0,  i  1 


\  all 


odd  values 


"'2  -  ^ 


gives  1  S  m^^  <  00  which  contributes  as 


pitl 

2p 


2  °°  [mj+(m^+2p)-2]/2 

‘(7)  1  2  '"l®!  ■  *  '"’'’2’ 


”1=1 

odd 


/4\*1  y  [-l)\  eo.  [m,(9, 

\txJ  2  Z/  "'1^1  ^  ^ 


2pe2] 


(3.14a) 


mj^=l 

odd 
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t 

i 


=  ■"  1 - 1  \  r  /  cos  [-m'(0 

\n/  2  zL  (-m|)[(-m^)  +  2p]  11 


0^)  -  2pe2] 


-m|=l 


odd 


=  ©'1  s 


m|=-l 


m;(ij^-'’2p)  ““  '"I*®! 


O2)  + 


(3.14b) 


odd 


Hence 


Sp  =  Sp  ^  "2p  "  Sp 


uu 

,  /1\  1  y  ,(-1)“^  ,  cos  [m,(6,  - 
\nf  2  Zw  nil  (nil  "  ^  ^ 


Og)  +  ^pBg] 


(3,15a) 


mi=-oo 


odd 


2  ,  /  ,  \P 
2  2p“ 


(1)^  i 


'cos  [nii(ei  -  62)  +  2pe2]  cos  [mi(Gi  -  Bg)  +  2p02]' 


nil  =-00 


nil  -  2p 


odd 


(3.15b) 


Further 


I 


cos  [(mi(Gi  - 


Oj)  Jpe^l 


m,  -  2p 


=  2 


cos  [(nil  -  2p)(ei  -  ©2)  +  2pGi] 


nil  -  2p 


nil  =-00 


nil  =-00 


odd 


odd 
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=  cos 


(2p(-l)  5! 


cos  -  2p)(Cj^ 


£'2 ) 


-  2p 


m^=-oo 

odd 


-  sin 


(2pfi^) 


sin  (m^  -  2p)((  ^ 


m^=-oo 


-  2p 


It  is  obvious  that 


s 


cos  nQ 
n 


=  cos  e(l  -  l)  +  cos  3( 


(I  -  i) 


n=-oo,  odd 


and ,  by  Ref .  28 , 


=  0 


I 

n=-oo 

odd 


sin  nO 
n 


=  2 


s 

n=l 

odd 


sin  nO 
n 


0  <  (j  <  n 


(3.16) 


Hence 


I 


cos[(m^(0j^  -  Og)  +  ZpOg] 


mj^=-oo 

odd 


-  2p 


-  2  2pGj^ 


0  <  C  <  n  (3.17; 
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Also,  proceeding  on  similar  lines, 


cos  -  Bg)  +  ^pe^] 


-  2 


0  <  G  <  n  (3.18) 


mj^=-oo 


Substituting  the  results  of  Eqs.  (3.17)  and  (3.18)  in  Eq.  (3.15b)  gives 


C„  =  - [sin  2p0„  -  sin  2pe,  ] 

2p  n  2p  2  I-* 


(-1)^  sin  p(e^  -  0^)  cos  p(0j^  f  Bg) 


0  <  (e^  -  Bg)  <  n 


(3.19a) 


p  ^  0  (3.19b) 


Equations  (3.8c),  (3.11),  and  (3.19b)  give 


r  1  A  ("1)^  P(®1  "  69)  6p) 

y(t)  y(t-i)  =  [i-f  (e,-e,)]  -  i  2  - 5-^2 - 1-^ 


0  <  (01  -  Bg)  < 


(3.20) 


D .  OUTPUT  PROCESS 


1.  Expansion  of  z(t, 


The  expansion  of  the  FM  pulse  train  can  be  rewritten  by  sub¬ 
stituting  the  result  of  Eq .  (3.20)  in  Eq.  (3.7),  which  gives 


.(t)  =7  (e,  -  ej)  +  2^  J 


(-1)^  sin  p(Ci  -  e^)  cos  p(ei  +  e 
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“  (-1)^  sin  p(Gj^  -  Og)  cos  p(Gj^  +  Bg) 

= «>i  -  02)  * "  Z - 5 

P=1 


=  7  "1  -  S)  +  I  <-*> 


sin  2pe,  -  sin  2pG„ 
p  _ i _ £ 


(3.21) 


where 


Gj^  =  total  phase  of  the  input  process  x(t)  at  time  t 

e  =  total  phase  of  the  input  process  x(t)  at  time  (t 
2 

It  may  be  noted  that 


e,  -  0„  =  oci  T  +  [\ir(t)  -  ilf(t  -  T)] 

X  A  C 


(3.22) 


01  +  02  =  2 


[“c  (^  '  I)  "  ^  ^ 


where 


\lr(t)  =  information  for  phase  modulation 


lir(t)  =  information  for  frequency  modulation 

It  may  be  added  that  all  the  expressions  derived  above  are  valid  for 
those  cases  where  x(t)  is  subjected  to  hard  limiting  with  wideband 
limiters . 

;; .  Autocorrelation  Function 

The  derivation  of  the  au toco irelat ion  function,  which  is  quite 
involved,  is  shown  in  Appendix  B.  The  final  result  is  contained  in  Eq. 
(B.19): 


R  (t,t'  ) 
z 


*  = (if  (“y 


“  ”  2/(2pP)  sin^  (2pu)  +  noj  )  |  1 

^  ^  _!l - - -  cos  (2paj^  +  '  ) 

P=1  n=-°°  ^  (B.19) 
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3 .  Power  Spectral  Density 


The  expressions  for  the  power  spectral  density  of  the  output 
process  are  developed  assuming  rectangular  pulse  trains,  arbitrarily 
shaped  pulse  trains,  and  power  spectral  density  with  no  modulation. 

a.  Rectangular  Pulses 

In  the  case  of  rectangular  pulse  trains  it  is  easily  seen 

that 


00  00  2  /  \ 

2J  (2pP) 


2,  Z  - 2 -  5j 


p=l  n=-oo 


8[f  -  (2pf  +  nf  )]  +  5[f  +  (2pf  +  nf  )] 
'■  '  c  m  '  c  m 


(3.24a) 


(1)^  .  (2, 


o  o  5(f  -  f  )  +  +  f  ) 

2  .  2  (juT\  ^  m^  '  m' 

2  /  2 


00  00  2/  \  2/  \T 

2J  (2pP)  sin  (2pu)  +  noj  )  — 

"  c  m  2 


S  I  ^ 

p=l  n=-oo 


&[f  -  (2pf^  +  nf^)]  +  6[f  +  (2pf^  +  nf^)] 


(3.24b) 


b.  Arbitrarily  Shaped  Pulses 

For  arbitrarily  shaped  pulse  trains  it  is  better  to  cast  the 
expression  for  into  a  more  transparent  form  as  follows; 


SEL-65-05G 


-  82  - 


t-rw  v.'v  1^*  «!>«  wMMit  MKHrMna  !•_■(  a.  m  n  ■  r  -v  m  r  »r  Jt.  'vj*  r_«  rjk  r_fl  n  m  -  rj  a*  »  "  ^  n  m  n 


9- 

i 


2  /  2  „  o  S(f  -  f  +  -  f 


i  S  i4^  )' 

^  2  222  '  c  m 

p=l  n=-oo  p  n  f  T 


6[f  -  (2pf^  +  nf^)]  +  S[f  +  (2p£^  +  nf^)] 


(3.25a) 


r  „  „  &(f  -  f )  +  6(f  +  f ) 

(4)^  <  (2f^T)^  6(f)  +  2t3^f^  (t  sine  fl)^  - - — ^ - — 


1  1 


2(2pf  +  nf  J^(2pP)  2 

- 5 - E - 2 -  (T  sine  fT)^ 


p=l  n=-oo 


&[f  -  (2pf^  +  nf^)]  +  &[f  +  (2pf^  +  nf^)] 


(3.25b) 


(  „  „  6(f  -  f  )  +  5(f  +  f  ) 

(4)2  (T  Sine  fT)2  /(2f^')2  6(f)  +  2p2f2  - EL_ - !L 


s  s 


2(2pf^  +  nf^)2  j2(2pp) 


p=l  n=-oo 


6[f  -  (2pf^  +  nf^)]  +  6[f  +  (2pf^  +  nf^)] 


(3.25e) 
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f  k  Fk  •  m  F_i«  n_k.  RJC  f\M  T-li  Ji  rj^  fVA  rj<  *J» 


•-F  ‘ji  -j«r>  ■'>  \aV» 


Also , 


jTl(t)  ->  T  sine  fT 


Hence 


B^(T) 


=  (if  Jrlit)  *  Jrlit)  * 


/  ^2  2  2 
(2f  )  +  2P  f  cos  o)  T 

c  m  m 


UO  00 


I 

p=l  n=-o 


2(2pf  +  nf  J^(2pp) 

'  c  m  n 


cos  (2pcD  +  ncD  )t 
c  m 


(3.26) 

Let  p(t)  be  any  arbitrarily  shaped  pulse  with  P(f)  as  its  transform, 
and  whose  width  is  such  that  no  overlapping  of  the  successive  pulses 
occurs.  Then  for  such  :•  pulse  train  S^(f)  becomes 


6(f  -  f  )  +  &(f  +  f 

S,(f)  =  (4)^  [P(f)]''‘  {{2tS  +  2p^f" - 


m 


00  00 


I  I 

p=l  n=-oi 


2(2pf  +  nf  J^(2pP) 

'  c  m  n' 


5[f  -  (2pf  +  nf  )]  +  6[f  +  (2pf  +  nf  ) ] 
^  '  c  m^-^  '  c  m' 


(3.27) 


which  gives  the  power  spectral  density  for  a  train  having  arbitrarily  but 
fixed  shaped  pulses  at  every  zero  crossing. 

c.  Output  Spectra  with  No  Modulation,  p  =  0 

It  can  easily  be  shown  from  Eq.  (3.25b)  that  the  spectral 
power  density  with  no  modulation  i.. 
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s^(l) 


no 

mod 


(4)^  (2f 


“  6(f  -  2pf^)  +  8(f  +  2pf^) 

5(f)  +  ^  2  sine  (2pf^T)  - 2 - 

p=l 

(3.28) 


The  above  result  can  easily  be  derived  in  a  simpler  way  by  noting  that 
the  output  with  no  modulation  consists  of  a  pulse  train  of  random  phase, 
the  repetition  rate  being  2f^  and  each  pulse  being  R(t). 

4.  Output  Voltage  Vs  Frequency  Characteristic 

It  is  convenient  to  express  the  sensitivity  of  the  system  shown 
in  Fig.  11  in  terms  of  the  dc  output  voltage  From  the  expression 

for  S  (f)  it  is  easy  to  show  that  the  rms  voltage  output  at  any  fre- 

Z 

quency  f^  is 


pv. 


dc 


f 


sine  f  T| 
m 


m 


(3.29) 


The  linearity  of  such  a  system  is  therefore  excellent. 

E .  INTERMODULATION  PRODUCTS 

From  the  power  spectral  density  of  z(t),  it  is  seen  that  the  inter¬ 
modulation  products  can  be  tagged  by  two  indices,  p  and  n.  Hence 
IM(p,n)  would  imply  the  following: 

1.  The  intermodulation  product  IM(p,n)  belongs  to  the  p  spectral 
zone,  p  assuming  all  positive  integer  values. 

2.  The  order  of  the  intermodulation  product  n  can  have  all  integer 
values  from  -<»  to  +<». 

3.  The  frequency  of  the  intermodulation  product  IM(p,n)  is 
(2pf  +  nf  ) . 

'em'  22 

4.  When  the  power  in  the  wanted  information  is  [2(3  sin  (aj^T/2)], 

then  the  power  in  the  intermodulation  product  IM(p,n)  is 

2J^(2pP)  sin^  (2pm  +  rei  )  | 

_ E _ ni — -  ( 3 . 30 ) 

2 

P 
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since  the  power  in  IM(p,n)  is  at  a  frequency  of 
is  possible  for  certain  values  of  p  and  negative  values  of  n  that 
IM(p,n)  will  lie  in  the  wanted  spectrum  and  hence  cannot  be  filtered 
out.  This  intermodulation  product  appears  in  the  output  as  an  inter¬ 
ference  at  a  frequency  which  is  not  harmonically  related  to  the  informa¬ 
tion  frequency  and  is  therefore  sometimes  referred  to  as  a  "spurious 
output."  In  such  cases,  it  is  more  convenient  and  meaningful  to  define 
a  ratio,  termed  the  "interference  ratio,"  or  IR(p,n). 

1 .  Interference  Ratio,  IR(p,n) 

The  interference  ratio  IR(p,n)  is  defined  as  the  power  in 
IM(p,n)  when  the  power  in  the  wanted  information  is  normalized  to  1. 
Hence 


IR(p,n)  =  interference  at  frequency  (2pf^  +  ’’^m^ 


(2pP)  sin  (2pcu^ 

oTt 

pP  sin  ~ 


+  nu) 

m 


2 


(3.31) 


which  can  also  be  expressed  in  decibels. 

Almost  all  existing  wideband  FM  systems  have  operating  parameters 

such  that  the  values  of  oj  T/2  and  (2pa)  +  nao  )(t/2)  are  extremely 

m  '  cm' 

small  and  the  sines  may  be  replaced  by  their  arguments.  This  simplifi¬ 
cation  gives 


IR(p,n) 


fj„(2pP) 
■p^ 


=  20 


log 


10 


pP 


db 


(3.32) 
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provided  f^^O,  P^O,  pe  (l,oo),  and  n  ^  0.  It  may  be  noted  that 
n  has  a  negative  integer  value  in  IR(p,n). 

The  output  z(t)  is  filtered  to  recover  the  original  information. 
Let  the  output  of  the  final  lowpass  filter  be  Hence  all  the  inter- 

modulation  components  which  satisfy  the  condition  (2pf^ 
are  filterable  and  are  therefore  of  no  consequence  as  far  as  the  overall 
system  performance  is  concerned. 

The  strongest  interference  usually  occurs  for  p  =  1 ,  and  for 
the  highest  negative  value  of  n  which  satisfies  the  condition 

(2f  +  nf  )  S  B. 

c  £p 

that  is,  n  is  the  Integer  which  satisfies 


n  ^ 


2f 


ip 


(3.33) 


Since  the  highest  negative  value  of  n  will  occur  for  the  highest 
in-band  frequency,  the  worst  interference  will  be  for  that  n  which  is 
obtained  when  f^  is  replaced  by  ^£p-  The  worst  interference  occurs  at 
a  frequency  (2f^  +  Has  a  ratio 


IR(l,n) 


=  20  log 


db 


(3.34) 


where 


=  modulation  index 


Usually  a  ±k  percent  system  Implies  that 


P  = 


100  f 


m 


0.01 


(3.35) 
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Therefore 


db 

(3.36) 

It  may  be  noted  further  that  only  the  first-order  interference  components 
need  be  computed,  since  in  most  of  the  wideband  FM  systems  the  total  con¬ 
tributions  due  to  higher  order  components  amount  to  less  than  10  percent 
of  the  first-order  term. 

2 .  Pulse  Shape  and  Interference  Ratio 

In  computing  the  first-order  effects,  the  interference  ratio  is 

calculated  at  a  point  where  the  equality  f  =  (2pf  +  nf  )  is  satisfied. 

m  cm 

Hence  the  interference  ratio  is  essentially  unaffected  by  the  shape  of 
the  pulses . 

F .  RESULTS 

1 .  Theoretical  and  Experimental  Results 

In  order  to  be  able  to  select  optimum  system  parameters,  given 
a  set  of  specifications,  Eqs.  (3.3l)  -  (3.36)  are  used  to  compute:  the 
interference  ratios  in  Table  8  and  the  curves  of  Figs.  12  and  13. 

1.  Table  8  shows  the  interference  ratio  in  decibels  for  various 

deviations.  The  information  bandwidth  is  assumed  to  be  0  to  1 , 
while  the  nominal-carrier  frequency  Is  set  at  2.  The  interference 
ratios  IR(1,-3),  IR(1,-5),  IR(1,-7),  and  IR(1,-11),  which 

cause  spurious  outputs  at  frequencies  1,  2/3,  l/2,  and  l/3,  are 
shown  in  the  table. 

2.  The  curves  of  Fig.  12  show  graphically  the  relationship  between  the 
interference  ratio  and  the  deviation  when  the  ratio  of  the  nominal- 
carrier  frequency  to  the  information  bandwidth  is  2. 

3.  The  curves  of  Fig.  13  show  the  behavior  of  spurious  output  vs 
carrier  frequency  at  constant  deviation 


IR(l,n)  =  20  log 


10 


IJ|  I [0.02  k  (f  /f  )] 

I  n  I  ^  c^  m^ 

0.01  k  (f  /f  ) 

'  m' 


ft) 
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TABLE  8.  INTERFERENCE  RATIOS  FOR  VARIOUS 
DEVIATIONS  AND  BANDWIDTHS 

Normalized  information  bandwidth,  BW  =  1 
Normalized  nominal  carrier  frequency,  f  /BW  =  2 
(All  values  in  minus  decibels) 


System 

Deviation 

(±k^) 

ir(i ,-ll) 

f/BW  =1/3 

IR(l,-7) 
f/BW  =1/2 

IR(1 ,-5) 
f/BW  =2/3 

IR(1,-3) 
f/BW  =  1 

0 

10 

— 

83.6 

43.6 

12.5 

— 

108.0 

— 

— 

20 

— 

— 

59.9 

31.8 

25 

— 

75.1 

52.4 

28.2 

30 

— 

— 

46.4 

25.2 

37.5 

— 

55.4 

— 

— 

40 

— 

— 

37.4 

20.9 

42 

77.1 

— 

— 

— 

50 

— 

42.4 

30.9 

17.8 

Note:  Values  not  indicated  were  either  negligibly  small  or 
not  computed. 


The  experimentally  measured  results,  which  are  also  shown  on 
Figs.  12  and  13,  are  in  close  agreement  with  the  computed  values.  For 
example,  spurious  output  measurements  which  were  made  for  various  devia¬ 
tions  and  for  nominal-carrier  and  input  frequencies  were,  in  every  case, 
within  2  db  of  the  theoretically  computed  values.  A  typical  example  of 
these  experimental  measurements  is  reproduced  below  from  data  obtained 
on  6  November  1964 : 


Conditions 


Nominal-carrier  frequency; 
Deviation : 

Test  gear: 


900  kc 
+40  percent 

> 

FM  multivibrator  modulator 
FM  pulse  train  demodulator^ 
Spectrum  analyzer 


prototypes 
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SYSTEM  DEVIATION  tk 

0  10  20  30  40  SO 


FIG.  12.  CURVES  SHOWING  THE  INTERFERENCE  RATIO  VS  DEVIATION 

FOR  FIXED  f  /f  =2. 

c'  nun 


NOMINAL- CARRIER  FREQUENCY  (normollztd) 


FIG.  13.  CURVES  SHOWING  SPURIOUS  OUTPUT  BEHAVIOR  VS  CARRIER  FREQUENCY 
AT  CONSTANT  DEVIATION, 
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In  Table  9  the  computed  and  measured  values  for  various  IR(p,n)  are 
compared . 


TABLE  9.  COMPARISON  OF  COMPUTED  AND  MEASURED  VALUES 
FOR  VARIOUS  IR(p,n) 


Input 

(kc) 

IR(p,n) 

Frequency 

(kc) 

Computed  Value 
(db) 

Measured  Value 
(db) 

500 

IR(1,-3) 

300 

-26.8 

-29 

IR(l,-4) 

200* 

1 

CO 

00 

-38 

400 

IR(l,-5) 

200* 

-52.5 

-53 

550 

-25.1 

-27 

250 

IR(1,-6) 

300 

-42.4 

-40 

IR(l,-7) 

50 

-71.3 

— 

IR(1,-14) 

400 

-69.8 

— 

100 

IR(l,-15) 

300 

-- 

— 

* 

The  input  frequency  had  to  be  rocked  slightly  to  distinguish 
between  these  two  components. 


The  linearity  and  dc  drift  measurements  were  made  on  a  typical 
FM  pulse  train  demodulator  with  a  nominal-carrier  frequency  of  900  kc  and 
a  deviation  of  ±40  percent.  The  test  gear  used  was  a  five-place  digital 
voltmeter,  a  frequency  counter,  and  a  PM  prototype  demodulator.  It  was 
found  that  the  maximum  deviation  from  best  straight  line  through  zero  was 
better  than  ±0.1  percent  and  that  the  dc  drift  was  better  than  0.01  per¬ 
cent  per  degree  Centigrade. 

2 .  Advantages  over  Conventional  Detectors 

As  far  as  low  carriers  and  wide-deviation  FM  systems  are  con¬ 
cerned,  this  demodulation  technique  affords  a  better  way  to  recover 
original  information  than  the  conventional  FM  discriminators  or  ratio 
detectors.  The  conventional  detectors  having  a  maximum  deviation  from 
linearity  of  less  than  ±1  percent  to  over  ±50  percent  of  the  center 
carrier  frequency  are  extremely  difficult  to  realize  [see  Ref.  24,  p. 

494;  Ref.  29,  p.  1104;  and  Ref.  30,  p.  214]. 
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1.  FM  pulse  train  demodulators  represent  at  least  an  order-of-magnitude 
improvement  in  the  linearity  of  frequency  vs  output  voltage. 

2.  This  technique  results  in  a  device  having  excellent  dc  drift 
characteristics . 

3.  It  is  cheaper  and  easy  to  implement. 

4.  Since  the  i-f  filtering  before  detection  leaves  lower  sidebands  of 
the  third  harmonic  of  the  carrier  frequency  (resulting  in  amplitude- 
modulation),  the  interference  is  too  severe  to  permit  satisfactory 
operation  of  the  wideband  FM  discriminator. 

5.  The  ..'_i*dpass  filter  must  be  quite  sharp  and  yet  must  have  excellent 
envelope  delay  characteristics,  which  make  the  filtering  after 
limiting  (prior  to  detection)  undesirable. 

These  advantages  have  resulted  in  the  production  of  FM  demodulators 
having  thk'  xollowlng  specifications; 


Maximum 


Center  Carrier 
Frequency  f^ 

Deviation 

{±1o  of  f^) 

Information 

Bandwidth 

DC  Drift 

(fo/^C) 

from  Linearity 

{±1o) 

1  to  1800  kc 
going  up  in 
octaves 

40 

dc  to  1  Me 
depending 
on  f 

c 

0.01 

0.1 

7  Me 

30 

dc  to  5  Me 

-- 

1 

15  Me 

30 

dc  to  10  Me 

— 

1 
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IV.  REDUCTION  OF  SPURIOUS  OUTPUTS  IN  LOW-CARRIER, 
WIPE-DEVIATION,  WIDE- BANDWIDTH  FM  SYSTEMS 


A.  FORMULATION  OF  THE  PROBLEM 

In  many  systems  and  practically  all  electromagnetic  recorders,  the 
channel  has  a  bandwidth  from  a  few  hundred  cycles  to  a  few  megacycles. 

In  most  cases  the  information  has  a  bandwidth  down  to  dc  and  therefore 
some  form  of  a  modulation  scheme  is  required.  The  bandwidth  of  the  chan¬ 
nel,  the  modulation  scheme,  the  tolerable  strengths  of  extraneous  outputs, 
and  the  signal- to-noise  ratio  are  some  of  the  factors  which  determine  the 
bandwidth  of  the  information  that  can  be  recorded.  Frequently  vestigial 
FM  recording  is  used  in  order  to  record  information  of  much  larger  band- 
widths  than  is  possible  with  conventional  FM.  The  spectrum  of  a  typical 
vestigial  FM  recording  is  shown  in  Fig.  14.  A  high  modulation  index  is 
needed  in  order  to  obtain  an  acceptable  signal- to-noise  ratio. 


MAXIMUM  EXCURSION  OF 


FIG.  14.  TYPICAL  BANDWIDTH  OF  CHANNEL,  USING  VESTIGIAL 
FREQUENCY  MODULATIONS. 

In  the  frequency-modulation  case,  some  relations  between  the  band¬ 
width  of  the  modulated  wave  and  the  frequency  deviation  are  given  by  Harris 
[Ref.  31]  and  Abramson  [Ref.  32].  It  is  assumed  in  the  analysis  by  Harris 
and  also  by  Abramson  that  the  input  is  a  narrowband  random  process.  This 
is  shown  in  Fig.  15. 

In  the  systems  analyzed  in  this  chapter,  the  narrowband  assumption  is 
not  valid.  A  typical  spectrum  of  a  wideband  system  is  shown  in  Fig.  17, 
p.  97.  The  spectrum  centered  at  f^  partially  overlaps  the  spectrum 
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g(t) 


FREQUENCY 

MODULATOR 


y(t) 


FIG.  15,  MODULATED  SPECTRUM  WHEN  »  f 

h  ymm 


centered  at  as  a  result  of  the  nominal -carrier  frequency  being 

comparable  to  the  information  bandwidth.  Because  of  this  partial  over¬ 
lapping,  "aliasing"  occurs  and  outputs  appear  in  the  demodulated  informa 
tion  which  cannot  be  attributed  to  a  harmonic  distortion  of  the  input 
information.^  These  so-called  "suprious  outputs"  vanish  under  certain 
obvious  conditions  which,  however,  are  unacceptable  for  the  reasons  indi 
cated  below: 


Condition 

1.  Sufficiently  high  nominal-carrier 
frequency  so  that  overlapping  of 
the  spectra  does  not  occur. 

2.  Reduced  information  bandwidth. 

3.  Reduced  deviation. 


Reason  for  Nonacceptance 

Because  of  transducer  electro¬ 
magnetic  limitations. 

Makes  the  channel  uneconomical 

Results  in  rapid  deterioration 
of  si gnal- to-noise  ratio. 


Wang  and  Wade  [Ref.  33]  discuss  this  phenomenon  in  connection  with 
superregenerative  parametric  amplifiers. 
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The  analysis  of  this  chapter  is  devoted  to  finding  acceptable  methods,  if 
there  by  any,  for  reducing  the  strengths  of  the  spurious  outputs. 


B.  ANALYSIS  OF  THE  PROBLEM 

It  is  assumed  that  the  information  g(t)  has  a  bandwidth  ("^^m’^mm^’ 

where  f  denotes  the  maximum  frequency  present  in  the  information, 
mm 

Define 


t(t)  = 


d| 


If  f,.  is  the  nominal-carrier  frequency  of  a  frequency-modulated  wave 
h 

v(t),  and  f,  »  f  ,  then  after  normalizing  the  amplitude  to  1, 
h  mm 


y(t)  =  cos  [ojj^t  +  +  ^(t)] 


(4.1) 


It  can  easily  be  shown  that  if  \lr(t)  is  a  gaussian  random  process  with 
a  finite  mean  square  value  [Ref.  32], 


Ry(T)  =  I  exp  [-R^(O)]  exp  [R^(t)]  cos  (4.2) 

where  R  (t)  is  the  autocorrelation  function  of  y(t).  The  power  spectral 

y 

density  of  y(t)  is  (see  Fig.  15), 


S  (f)  =  ^  exp  [-R^(O)]  n{exp  [R^(t)]}  *  [5(f  -  f^^)  +  6(f  +  f^^)] 

(4.3) 


where  the  Fourier  transform  of 


is 


n(exp  [R^(t)]]  =  J  exp  [R^('r)]  exp  (-icux)  dT 


Although  the  information  g(t)  is  bandlimited  to  f  l^mm^’ 

seen  from  Fig.  15  that  the  bandwidth  of  ®y(^)  given  by  expression 

(4.3).  In  the  following  discussion  it  is  assumed  that  ®y(f)  exists  in 
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a  band  2f  centered  at  ±f,  and  is  practically  zero  elsewhere, 
ymm  h 

Usually  >  f__. 

ymm  mm 

If  the  conditions  in  Fig.  15  exist,  that  is,  then 

demodulation  by  any  standard  technique  will  result  in  the  recovery  of 

the  information  g(t)  without  spurious  outputs. 

If  y(t)  is  now  translated  down  to  another  frequency,  say  fj^  is 

greater  than  f  as  shown  in  Fig.  16,  the  original  message  can  still 

®  ymm 

be  recovered  without  any  spurious  outputs.  This  is  evident  from  the 
fact  that  if  one  reverts  back  to  the  original  carrier  frequency  fj^  by 
translation,  as  shown  in  Fig.  16,  one  obtains  y(t)  exactly.  Hence  th6 
whole  process  is  "reversible." 


y(t) 

MULTIPLIER 

BANDPASS  FILTER 
BAND  CENTER 

i(t) 

> 

( 

BANDWIDTH 

CO»[(W^-W^.)f] 

z(t) 

MULTIPLIER 

BANDPASS  FILTER 
BAND  CENTER 
BANDWIDTH 

ylt) 

> 

( 

CO*[(W^-W^.)t] 

FIG.  16.  MODULATED  SPECTRUM  WHEN  f  ’  §  f 

h  ymm 
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The  case  where  y(t)  is  translated  down  to  a  carrier  frequency 

less  than  f  is  shown  in  Fig,  17.  This  condition  is  usually  pre- 

ymm 

valent  in  wideband  FM  recorders.  Partial  overlapping  of  the  positive  and 

negative  spectra  occurs  and  spurious  outputs  result.  The  shape  of  the 

spectrum  of  S  (f)  is  different  from  S  (f);  and  given  such  a  z(t), 
z  y 

one  cannot  revert  back  to  y(t)  by  translation.  The  process  is  there¬ 
fore  "irreversible." 


C.  INTERRELATIONSHIPS  OF  f  ,  f  ,  and  f  FOR  FILTERABLE 

_  c  ymm  Ip 

SPURIOUS  OUTPUTS  ^ 

If  one  demodu.I  ates  z(t)  and  passes  the  output  through  a  lowpass 
filter  to  recover  the  message  g(t),  a  portion  of  the  spurious  outputs 


y(») 

MULTIPLIER 

LOWPASS  FILTER 

2(t) 

X 

cos 


FIG.  17.  MODULATED  SPECTRUM  WHEN  f. 


C  fymm  (SHOWS 
y(t)  cannot 


ALIASING).  As  a  result  of  aliasing, 
be  recovered  from  z(t).  The  dotted  lines  in  the 
lower  figure  indicate  the  shape  of  S  (f)  and  all 


,(f)  for  f  >  f 


ymm 
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generated  due  to  the  aliasing  phenomenon  is  filterable  (Fig.  18).  To  be 
more  specifi 
Fig.  17  and 


more  specific,  if  the  lowpass  filter  has  a  bandwidth  of  f ,,  ,  then  from 

ip 


i  [f  -  (f  -  f  "I]  =  2f  -  f 
ip  c  ymm  c  c  ymm 


(4.4) 


minimum  distortion  occurs  in  g(t)  since  the  spurious  outputs  are 

filterable.  Hence  if  g(t)  is  bandlimited  to  frequencies  less  than 

(2f  -  f  ),  good  demodulation  is  possible. 

'  c  ymm' 


1  AU/DAQQ  Pit  TPD  1 

g(t) 

FIG,  18.  LOWPASS  FILTERED  OUTPUT. 

In  the  frequency-modulation  case,  the  spectrum  theoretically  has  an 

infinite  bandwidth.  In  all  practical  systems,  however,  S  (f)  is  sig- 

z 

nificant  over  a  small  frequency  interval.  All  components  of  S^(f)  out¬ 
side  this  interval  may  be  vary  much  below  the  signal- to-noise  ratio  of 
the  system,  which  sets  a  criterion  in  estimating  from  Eq .  (4.4). 

To  illustrate  the  application  of  the  result,  consider  a  simple  case  of 
monotone  modulation.  If  the  modulation  index  P  for  the  maximum  frequency 

of  information  f  is  such  that  only  p  lower  sidebands  are  signifi- 
mm 

cantly  different  from  zero,  then  for  the  case  of  no  spurious  outputs  in 
the  filtered  output  due  to  aliasing, 

f  S  (2f  -  pf  ),  that  is,  f  =  { — — r)  f  (4.5) 

mm  c  mm''  mm  \p  +  1/  c 

D.  REDUCTION  OF  SPURIOUS  OUTPUTS 

Since  the  bandwidth  of  the  channel  is  less  than  f  ,  the  best  one 

ymm 

can  do  for  minimization  of  spurious  components  is  to  make  ^^(f)  in 
Fig.  17  coincide  with  ®y(^)  throughout  the  entire  bandwidth  of  the 
channel,  as  shown  by  the  dotted  lines  in  Fig.  17.  This  approximation 
can  be  done  by  appropriate  filtering  or  phasing  techniques. 
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1.  Filtering  Technique 


It  is  clear  that  the  filtering  should  be  performed  when  the 
conditions  in  Fig.  15  prevail.  The  bandpass  filter  must  have  an  approxi¬ 
mately  constant  group  delay  characteristic  close  to  the  first  null.  The 

slope  of  the  filter  depends  upon  f  ,  f  ,  f .  ,  the  modulation  index 
^  mm  c  n 

P,  the  signal-to-noise  ratio  of  the  system,  and  the  rate  of  cutoff  of 
the  final  lowpass  filter  at  the  output  of  the  receiver.  With  a  set  of 
Bessel  function  tables  and  a  reasonably  low  value  of  f^^,  the  system 
parameters  can  be  chosen  so  as  to  result  in  filters  that  are  realistic 
and  inexpensive.  The  complete  system  which  implements  the  filtering 
technique  is  shown  in  Figs.  19  and  20. 

2.  Phasing  Technique 

The  block  diagram  which  implements  the  phasing  technique  is  shown 
in  Fig.  21.  The  modulator  of  Fig.  21a  and  the  demodulator  of  Fig.  21b 
have  both  resulted  from  the  following  theoretical  analysis: 

y(t)  =  cos 

=  cos  (u)^t  +  [cos  t(t)]  -  sin  (oD^t  +  0^)  [sin  t(t)] 

Taking  Fourier  transforms , 

Y(f)  = 

i<ti  -itti 

c  c 

-  6(f  -  f  )  -  ^  6(f  +  f  )  *  rUsin  t(t)} 

C  « 1  c 

J  (4-6) 

Define 

C  .  =  :Hcos  il/(t)  j  *  S(f  +  f^)  (4.7a) 

-  X  c 

c 


i  <ti 


—  8(f  -  f,)  t 


-i(}j 


*  <'l  [cos  '1^(  t )  ] 
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FIG.  19.  BLOCK  DIAGRAM  AND  SPECTRA  OF  THE  PROPOSED  SYSTEM 
USING  THE  FILTERING  TECHNIQUE. 


cos  ZTTfjt 
(f^  ARBITRARY) 

FIG.  20.  DEMODULATION  PROCESS  USING  THE  FILTERING  TECHNIQUE. 
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a.  Modulation 


b.  Demodulation 

FIG.  21.  BLOCK  DIAGRAMS  OF  THE  PHASING  TECHNIQUE. 


C^^  =  J3{cos  t(t)J  *  6(f  -  f^) 

c 

=  3{sin  *  6(f  +  f^) 

c 

S^^  =  :3(sin  ^If(t)J  *  &(f  -  f^) 

C 

[F(f)]^  =  F(f)  I  [I  +  sgn  f] 

=  F(f)  for  all  f  <  0 

=  0  for  all  f  >  0 


(4.7b) 

(4.7c) 

(4,7d) 


(4.8) 
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[F(f)]_  ^  F(f)  i  [1  -  sgn  f] 


=  F(f)  for  all  f  <  0 
=  0  for  all  f  >  0 


Hence 


Y(f) 


i<$i  -iij)  i'^  lit 


also  define 


r(t)  =  sin  +  41^  +  \lr(t)] 


R(f)  = 


14  -14 


*  3{cos  'lf(t)} 


1 4  -14 

■  f^)  .  £^6(f  *  t^) 


■X  3(sin  \lA(t)j 


14  -14  14  -14 

c  c  c  c 


Defining  the  Hilbert  transform  of  x(t)  as  x(t), 


x(t)  =  Cauchy's  principal  value  of 


JOO  , 

t  - 

IJO 


ds 


(4.9) 


(4.10) 


(4.11) 


(4.12) 
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»> 


(4.13) 


Using  the  relation 


F(f)  =  [F(f)]^  +  [F(f)] 


(4.14) 


and  subtracting  R(f) 


from  Y(f)  gives 


Y(f)  -  R(f) 
2 


(4.15) 

A  comparison  of  the  above  equation  with  Eq.  (4.10)  for  Y(f) 
shows  that  the  wanted  output  is  precisely  Y(f)  -  R(f)  since  the  inter¬ 
fering  portions  of  the  positive  and  negative  spectra  nullified  each  other. 
The  block  diagram  of  the  modulator  based  on  this  analysis  is  shown  in 
Fig.  21a.  To  recover  the  message,  the  systems  shown  in  Fig.  19  or  Fig. 
21b  can  be  used. 

It  may  be  added  that  Hilbert  transformers  can  be  approximated  to 
a  desired  accuracy  over  a  wide  frequency  interval  by  several  methods, 
e.g.  ,  an  active  network,  a  passive  network,  or  a  tapped-delay-line  reali¬ 
zation.  Detailed  explanations  of  these  methods  are  contained  in  Refs.  34 
through  39.  The  phasing  method  with  modifications  can  be  used  to  develop 
complete  systems  using  SSB-FM,  SSB-PM  or  SSB-Hybrid  modulation.  Some 
experimental  results  on  SSB-FM  are  contained  in  Ref.  39. 

The  analysis  presented  in  this  chapter  has  resulted  in  techniques 
which  suppress  the  spurious  outputs  that  result  from  the  interaction  of 
the  positive  and  negative  spectra  when  the  nominal-carrier  frequency  is 
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sufficiently  lowered.  Block  diagrams  of  the  proposed  systems  are  pre¬ 
sented  in  Figs.  20,  21a,  and  21b.  Although  the  solutions  are  equally 
effective,  for  reasons  of  economy  the  two  systems  shown  in  Figs.  20  and 
21a  would  be  preferable. 
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V.  CONCLUDING  DISCUSSION 


An  analysis,  together  with  experimental  verification  of  the  theo¬ 
retical  results,  has  been  presented  which  allows  the  computation  of  the 
strength  in  any  output  component  and  the  "distortion"  (i.e.,  the  sum  of 
the  power  in  all  the  intermodulation  components  in  the  first-harmonic 
band)  when  many  signals  in  noise  are  subjected  to  hard  limiting.  Typical 
situations  in  which  these  results  have  proved  to  be  most  useful  are  as 
follows; 

1.  When  several  modulated  or  unmodulated  carriers  are  subjected  to 
hard  limiting--a  problem  of  common  occurrence  in  a  satellite  com¬ 
munication  systeiii  or  in  radar  problems  involving  multiple  echoes. 

2.  When  a  frequency-modulated  signal,  after  being  passed  through  a 
channel  having  arbitrary  output  vs  frequency  characteristics,  is 
subjected  to  hard  limiting. 

3.  This  study  may  be  used  by  various  organization,  such  as  NASA's 
Inter-Range  Instrumentation  Group  (IRIG),  for  drawing  up  specifi¬ 
cations  for  various  frequency-modulated  instrumentation  systems. 

For  the  general  case,  expressions  are  developed  which  give  the  signal- 
suppression  effect,  the  output  signal  power  distribution,  and  the  power 
in  the  various  crossproducts  due  to  hard  limiting.  A  few  results  that 
are  well  known  for  the  cases  of  one  signal  in  noise  and  two  signals  in 
noise  are  shown  to  be  valid  for  the  general  case. 

It  is  also  shown  that  the  behavior  on  hard  limiting  of  the  many- 
signals-in-noise  case  approaches  that  of  the  one-strong-signal- in-noise 
case  if  the  strong  signal  power  is  much  greater  than  the  sum  of  the  power 
in  all  the  remaining  weak  signals  and  noise.  It  should  be  stressed  that 
the  power  in  each  of  the  weak  signals  or  noise  need  not  be  equal  for  the 
above  approximation  to  be  valid.  Further,  the  case  of  hard  limiting  of 
many  signals  in  noise  is  reducible  to  that  of  two  signals  in  noise  if  the 
power  in  each  of  the  two  signals  is  much  greater  than  the  sum  of  the  power 
in  the  remaining  weak  signals  and  noise.  The  approximation  is  valid  for 
any  arbitrary  weak  signal  and  noise  power  distribution  since  it  is  the 


105 


SEL-65-056 


sum  of  the  power  in  them  that  is  of  consequence.  Also  the  case  of  p 
signals  in  noise  can  be  approximated  as  the  case  of  (p  +  l)  signals 
without  noise,  as  far  as  the  behavior  due  to  limiting  is  concerned. 
Further,  the  output  power  in  the  various  components  of  the  limited  signal 
is  relatively  independent  of  the  input  strong-signal-to-noise  ratios,  for 
ratios  larger  than  10  db.  For  the  case  of  p  equal  signals  (with  and 
without  noise),  the  results  of  hard  limiting  for  a  large  range  of  p  are 
presented. 

Some  of  the  contributions  of  this  research  are  listed  below: 

1.  A  means  of  predicting  the  approximate  behavior  on  hard  limiting  of 
many  signals  and  noise  has  been  developed. 

2.  Approximate  expressions  for  the  power  in  the  various  output  signal 
and  rrossproduct  components  for  any  p,  derived  analytically  are 
given.  These  expressions  give  results  which  are  very  close  to 
those  obtained  through  extensive  digital  computer  computations, 
the  error  in  most  cases  being  negligible. 

3.  The  judicious  selection  is  permitted  of  frequency  spacings  of 
several  modulated  carriers  so  as  to  minimize  the  distortion  over 
any  frequency  band  due  to  hard  limiting. 

4.  The  strength  of  the  interfering  components  present  at  the  output 
of  an  FM  system,  accomplishing  demodulation  by  hard  limiting  fol¬ 
lowed  by  integrating  an  FM  pulse  train  with  arbitrary  pulses  at 
every  p^”  crossover,  can  be  easily  evaluated  by  the  results  con¬ 
tained  in  this  report . 

5.  The  phenomenon  of  "aliasing,"  is  analyzed,  the  optimum  selection 
of  operating  parameters  is  permitted,  and  systems  are  proposed 
which  minimize  such  distortion. 
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APPENDIX  A.  USEFUL  INTEGRALS  OF  BESSEL  FUNCTIONS 


Recently  much  work  has  been  done  in  the  evaluation  of  integrals 
involving  Bessel  functions  and  products  of  Bessel  functions.  Some  of 
the  important  references  and  pertinent  equations  are  given  below. 


(-p^t^)  t^'^  J^(at)  dt 


(a/2p)-  P  (^) 

2p‘^  r  (v  +  1) 


i"i 


4  +  V 


V  +  1; 


exp  (-a^/4p)  (a/2p)^  T  (  ) 

2p^  r  (v  +  1) 


with  the  conditions 

R(n  +  v)  >  0, 


V  +  1-, 


(A.l) 


R(p^)  >  0 


For  tables  of  the  confluent  hypergeometric  function,  ^F^ ,  refer  to 
Slater  [Ref.  12]. 


exp  (-p^t^)  J^(at)  J^(bt)  dt 


I  (a/2p)>/2p)''  r  ^ 

r  (m  ♦  1)  r  {v  t  1) 


•I 


(-)^a/2p)2‘^  ( 


|.x  +  V  +  A  + 


k:(ii  +  1), 
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2Fi 


(A. 2) 


with  the  conditions 


r(h  +  +  a)  >  -1,  r(p  )  >  0 


where  generalized  factorial  function  defined  as  [Ref.  13] 


(«)„  =  n  (a  +  k  -  1) 


=  n.{a  +  1  )  (a  +  2 )  .  .  .  (a  +  n  -  1 )  n  g  l  (A .  3 ) 


(a)^  =1  a  ^  0 


The  following  Weber-Schaf heitlin  type  integrals  are  treated  in  great 
detail  in  many  references  [see  for  example  Watson,  Ref.  11] 


t  J  (at)  J  (bt)  dt  = - j 

•'o  ^  2r  (v  +  1)  r  ( 


(b/a)''(a/2)''-l  r  (H-LJl-lAiJ:) 


-  V  -  K 
2 


/  4  +  V  -  A  +  1  V  -  4  -  A+  1  .  ,  ^ 

2  1  \  2  ’2  ’  ’2 
\  a 


(a . 4a) 


with  the  conditions 


R(|i  +  V  -  A  +  1)  >  0,  R(A)  >  -1,  0  <  b  <  a 
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For  integrals  involving  the  product  of  more  than  two  Bessel  functions, 
it  is  possible  to  expand  the  other  Bessel  functions  in  powers  of  t,  and 
to  use  expressions  (A. 4b)  or  (A.l).  It  was  found  best  not  to  use  the 
above  expressions  for  the  general  case  as  it  is  easier  and  more  convenient 
to  carry  out  the  integrations  numerically  on  a  digital  computer. 
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APPENDIX  B.  DERIVATION  OF  AUTOCORRELATION 
FUNCTION  OF  THE  OUTPUT  z(t) 


It  Is  convenient  to  define 


5^(t)  =  [t(t)  -  i(t  -  t)] 


(B.l) 


Therefore,  (6^^  -  B^)  =  The  autocorrelation  function  of  the 

output  z(t)  is  defined  as 


=(i)  E 


V  *  *  2 

P=1 


(sin  2pe^^ 


sin  ZpG^^) 


Oj 


T  +  l^(t')  +  ^  (sin  2qe^^,  -  sin 


q=l 


(B.2a) 


where  E  denotes  the  expectation  and  0^^  denotes  the  value  of  6^^  at 
the  Instant  t.  Hence 


V*’'  •  •  V*  ’’’  *  4 i  i  ••••  ■  “  ••••’ 

■ 

.  .,t  ilytl  .  yi  l  .  j  IM.  ~  •».! 

.  j  tli  I—  •  «.  w*.!! 

•  *{1  V  >“  •  “  ■>•••4  ”{i  ^  ’  ■  •“  *•••}  ( B .  2b  ) 
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Now, 


E[|^(t)]  =  E[^it)  -  t(t  -  T)]  =  0 
E[|^(t')]  =  0 

E(sii'>  2p0^^)  =  E^sin  {2p[ui^t  +  \lf(t)]  +  2p(t>^}^ 

For  <|i  unifomly  distributed  on  [0,2jt],  and  il/(t)  and  it>  independent, 
c  c 

E[sin  2p0  ]  =  E(sin  2p[cD  t  +  'lf(t)]}  E[cos  2p(ti  ] 

X  w  c  c 

+  E{cos  2p[cD  t  +  '^(t)])  E[sin  2p(t>  ] 
c  c 

=  0 


Similarly, 

E(sin  2p02^) 

E(sin  2p02^, 

E(sin  2p0^^, 

Furthermore,  by  the  independence  of  the  random  process  ilf(t)  and  the 
random  variable  <t> 

c 

E[|^(t)  sin  2pe^^,]  =  E|t(t)  sin  {2p[a)^t'  +  +  2p<l'^]j 

-  EftCt  -  T)  sin  {2p[a)^t'  +  \lf(t')]  +  2p(t>^} 


=  0 
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which  further  results  in 


E{lT,(t')  [sin  -  sin  2pe5,^l) 


2t- 


>  =  0 


E{^^(t)  [sin  2pe^^,  -  sin  2pe^^, ]} 


2t' 


J 


Therefore , 


=  (l) 


+  E 


00  00 


S  I  ^ 

Lp=1  q=l 


•  (sin  2qej^^,  -  sin  2qe2^,) 


(B.3) 


Further  use  may  be  made  of  the  fact  that 


E[sin  2pej^^  sin  2qe^^,]  =  ^  E[cos  (2pe^^  -  2qe^^,)] 


It 


It  ^  It' 


-  I  E[cos  (2pej^^  +  2qe^^,)] 


=  E  cos  {2pa;  t  +  2p0  +  2pt(t) 

2  L  ^  ^ 

-  [2qa)  t'  +  2q4i  +  2q\l/(t')]} 

C  C  J 


-  E(cos  [2pa)  t  +  2p<)>  +  2p\l;(t) 

2  C  C 


+  2q(ju  t'  +  2q(t  +  2qilf(t')]] 
c  c 


=0  if  p  ^  q 
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Also,  for  all  Integral  value  of  p  and  q  e  (l,“) 


I  E[cos  +  2qe^^,)]  =  0 


The  other  terms  can  be  evaluated  in  a  similar  manner.  Hence, 


R^(t,t') 


(t,t')  +  E  < i  ^  [cos  2p(e^^  - 


p=l 


+  cos  2p(e2^  -  021'^  '  ”  ®lt'^ 


cos  2p(ej^^  -  02^,  )], 


(B.4) 


Define  t  -  t'  -  T,  that  is,  t'  =  t  -  T 


E[cos  2p(e^^  -  Qi^i)]  =  E{cos  2p[cD^t  +  \lr(t)  -  oo^t'  -  \lf(t')]} 


=  e|^cos  2p{cD^T  +  [\|/(t)  -  \l;(t  -  t)]}J 

=  E{cos  2p[(J0  T  +  ^_(t)]}  by  definition  of  I  (t) 

C  TT  I* 

=  EjRe  exp  (i[2pa)^T  +  2p|^(t)])j 

=  E^Re  {exp  [i2pco^T]  exp  [i2p|^(t)]}j 

=  Re  (exp  (i2pcD  t)  E[exp  (i2pt  (t)]]  (b.5) 

C 

By  definition  of  the  characteristic  function, 

'*’2p5  (f)  =  E[exp  (iaj2p5^)] 
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Here  is  is  assumed  that  \^(t)  is  stationary.  Therefore  =  'l^(t)  - 

ilr(t  -  t)  Is  also  stationary  and  depends  on  T  only.  Hence 


*2plJu)  ' 


(B.6) 


where  ,  is  the  characteristic  function  of  2p|  .  Therefore 
2p^,  T 

E[cos  2p(e^^  -  e^^,)]  =  Rejlxp  (i2pa.^T)  «>2p^(i5^)) 


(B.7) 


Similarly, 


E[cos  2p(e^^  -  Bgtt)]  =  E[cos  (2paj^T  +  2p^(t  -  T)}]  (B.8a) 


E[cos  2p(e2^  -  =  E[cos  {2poo^(T  -  T)  +  2p|^_^(t  -  T)}]  (B.8b) 


T-T' 


E[cos  2p(e^^  -  Bg^,)]  =  E[cos  (2paj^(T  +  T)  +  2p^^^^(t)}]  (B.8c) 


This  technique  allows  computation  of  any  specific  term,  given  t(t). 

To  proceed  further  a  specific  case  is  assumed  where  ^IfCt)  is  a 
random  process: 


\lr(t)  =  P  sin  (oD  t  +  0  ) 
'  m  m' 


(B.9) 


where  <f>  is  a  random  variable  with 
m 


(■t  '  to, 2.] 


p(»„)  =  { 


^^0  elsewhere 
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The  following  expansions  by  Jacobi  and  Anger  are  useful  [Watson, 

Ref.  11,  p.  22]: 

00 

cos  (z  sin  0)  =  Jq(z)  2  ^  (B.lOa) 

n=l 

00 

sin  (z  sin  0)  =  2  ^  '^2n+l^^^  ^  (B.lOb) 

n=0 

00 

cos  (z  cos  0)  =  Jq(z)  +  2  I  (  l)  (B.IOc) 

n=l 

00 

sin  (z  cos  0)  =  2  ^  (■!)"  '^2n+l^^^  ^  (B.lOd) 

n=0 

Also,  extensive  use  is  made  of  the  fact  that  for  0  uniformly  distributed 
on  (0,2n)  and  with  f(t)  and  0  independent  [Davenport  and  Root,  Ref.  7] 

E{cos  [f(t)  +  0]}  =  0 

Hence , 

E[cos  2p(0^^  -  e^^,)]  =  eJcos  {2p[a)^T  +  l.^(t)]3j 

=  cos  2pCD^T  E[cos  2p^^(t)] 

-  sin  2pa)^T  E[sin  2p|^(t)]  (B.II) 
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Now, 


E[cos  2p6  (t)]  =  E{cos  [2pP  sin  (ui  t  +  ♦  )  -  2pP  sin  (co  t  +  $  -  oi  t)]] 

T  mm  m  m  m 

=  Efcos  [2pP  sin  (u)  t  +  <ti  )]  cos  [2pP  sin  (ou  t  +  lO  -  t)]] 

m  m  m  m  m 

+  Efsin  [2pP  sin  (u)  t  +  ♦  )] 

m  m 

•  sin  [2pp  sin  (co  t  +  $  -  to  t)]} 

m  mm 

Using  the  above  expansions  and  taking  expectations, 


E[cos  2p5^(t)] 


=  E 


Jq(2pP)  +  2 


}  J-  (2pP)  cos  2n(a)  t  +  <t>  ) 
Zrf  2n  m  m 

n=l 


Jo(2pe)  t  2  ^  J2„(2pf5) 
m=l 


•  cos  2m(a)  t  +  4> 
m  m 


+  E 


2 


J„  ,(2pP)  sin  (2n  +  l)(ai  t  +  ((I  ) 
2n+l  '  m  m 


00 

2  7  J„  _(2pP)  sin  ( 2m  +  1 )  (oj  t  +  $ 

2m+l  m  m 

m=0 
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2  V 

Jq  (2pP)  +  0+0  +  4  y  - ^ -  cos  Snuj^T 

n=l 


^  jL+i(2pP) 

+  47  - - -  cos  (2m  +  l)(cD  t) 

/ .  2  '  m 

m=0 


=  Jq(2pP)  +  2  2  J^(2pP)  cos 


=  2 


n=l 


J  (2pP)  cos  no)  T 
n  m 


n=-oo 


(B.12) 


Also , 


E[sin  2p^^(t)]  =  E{sin  [2pP  sin  (oj^t  +  (t>^)  -  2pP  sin  (oD^t  + 


=  Efsin  [2p6  sin  (od  t  +  0  )]  cos  [2pP  sin  (cjd  t  +  <t  -  oj  t)]] 


-  Efcos  [2p6  sin  (ou  t  +  <1)  )]  sin  [2pP  sin  (cjut+it  -cot)]] 
'■  n,'-'  '  ni  mm 


=  E< 


2  /  J„  ,(2pP)  sin  (2n  +  l)(aLi  t  +  t 

Z-i  2n+l  m  r 


n=0 


J-(2pP)  +  2  7  J„  (2pP)  cos  2n(co  t  +  it>  -  O)  t) 

ov  I'K/  ^  2m'  '  '  m  m  m  ' 


m=l 
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-  E< 


Jo(2pP)  +  2  ^  J2|jj(2pP)  cos  2n(a)^t  +  <t^) 


m=l 


2  /  Jo„,  1  (2pP)  sin(2n+l)(ajt+(t  -cox) 
^  2n+l  m  m  m 


n=0 


=  0 


(B.13) 


Substituting  (B.12)  and  (B.13)  in  (B.ll)  gives 


E  cos  (2p[cD  T  +  £  (t)]}  =  y  J^(2pp)  cos  noj  x  cos  2p6D  X 
c  T  n  me 


n=-oo 


=  2 


J  (2pp)  cos  (2pco  +  nco  )x 


(B.14) 


n=-oo 


It  is  obvious  that 


E[cos  2p(e2^  -  02t'^^  "  E{cos  [2pcJD^X  +  2p^^(t  -  T)]} 


E{cos  [2p(jD^X  +  2p|_^(t)]} 


=  2 


J  (2pP)  cos  (2pcu  +  ncu  )x  (B.15) 
n  c  m  '  ^ 


n=-oo 
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E[cos  2p(e2^  -  6^)]  =  E{cos  [2paj^(T  -  T)  +  2p^_^(t  -  T)]} 


=  cos  2pa)^(T  -  t)  E[cos  2p|^  ^(t  -  T)] 

-  sin  2pa)^(T  -  T)  E[sin  2p|^_^(t  -  T)] 
It  is  evident  from  the  above  that 


E[cos  2p(e2^  -  )]  =  I  Jn^^PP)  cos  (2pco^  +  na)^)(T  t) 

n=-oo 

(B.16) 

Similarly,  it  can  be  shown  that 


E[cos  2p(0^^  -  02^,)]  =  I  J^(2pP)  cos  (2pa)^  +  n(jD^)(T  +  t) 

n=-oo 

(B.17) 

Hence 

E{cos  2p(0j^^  -  0^^,)  +  cos  2p(e2^  -  ©gt'^  '  ^^^®2t  '  ®lt'^ 


-  cos  2p(0j^^  -  02^.1  )]  =  ^  J^(2pP)[2  cos  (2paj^  +  nco^)T 

n=-oo 

-  cos  (2pa)  +  nCD  )(t  -  t)  -  cos  (2p(jL>  +  ncu  )(t  +  x)] 

'  *  c  m  '  c  m^' 


00 

=  2  ^  J^(2pP)[l  -  cos  (2pa3^  +  naj^)T] 

n=-oo 

•  cos  (2pa)  +  noD  )t 

'  c  m 
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00 


=  2 


y  J^(2pp)  2  sin^  (2paj  +  noj  )  | 
n  cm2 


n=-oo 


•  cos  (2pa)  +  ncjo  )t 

c  m' 


Also , 

R  (t.f)  =  E{|^(t) 

5t 

=  E{[ilf(t)  -  t(t  -  T)]  [\lf(t')  -  -  T)]} 

=  E[t(t)  -  E[i|r(t)  il/(t'  -  T)] 

-  E[\Jr(t  -  t)  t(t')]  +  E[\lr(t  -  T)  \|;(t'  -  T)] 
=  2R^(t)  -  [R^(t  +  T)  +  R^(t  -  T)] 

Taking  the  Fourier  transform  gives 


S^^(f)  =  2S^(f)  -  S^(f) 

=  2S^(f)  [1  -  cos  oiT] 

=4S^(f)sin^^  (B.18) 


Hence , 


R  (t.f)  =  R  (t) 

z  z 


I 


n=-oo 


2/(2pP) 


2 

P 


2 

•  sin 


+  no) 

m 


cos  ( 2paj 

c 


no)  )t 

m 
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Since  =  P  sin  (oo^t  +  it  can  be  concluded  that 


R^(t)  =  E[p  sin  (o.^t  +  0^)  P  sin  [o^^^t  +  -  o^^t)] 


=  ^  cos  OJ  T 
2  m 


which  implies 


„2  6(f  -  f  +  &(f  +  f 

S  -  g-  — _ 01- - ^ - 0!- 


Hence  in  this  case,  autocorrelation  of  the  output  process, 


p=l  n=-oo 


•  cos  (2pa)  +  no)  )t 

when  the  input  random  process  'lf(t)  =  P  sin  (oj^t  + 

Eq.  (B.19)  is  the  result  used  extensively  in  Chapter  III 

report . 


-  T)] 


(B.19) 


of  this 
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Attn:  SEPIR 


Chief  of  Naval  Research 
Navy  Dept . 

Washington  25,  D.C. 

2  Attn;  Code  427 
1  Code  4  20 

1  Code  463 

U.S.  Army  Electr.  Labs.  Comm. 
Mt .  View  Office 
P.  0.  Box  205 
1  Mt.  View,  Calif. 

Commanding  Officer 
ONR  Branch  Office 
1000  Geary  St . 

J  San  Francisco  9,  Calif. 

Chief  Scientist 
ONR  Branch  Office 
1030  E.  Green  St. 

1  Pasadena,  Calif. 


U.S.  Naval  Weapons  Lob. 
Dahlgren,  Va . 

1  Attn:  Tech.  Library 

Naval  Ordnan< e  Lab. 

Corona .  California 
I  Attn:  Library 
I  H.  H.  Wieder,  423 

Commanding  Officer  (ADL) 
USN  Air  Dev.  Ctr. 

1  Johnsvtllc.  Pa.  18974 

Commander 

USN  Missile  Center 
Pt,  Mugu,  Calif. 

1  Attn:  N03022 


Commandant 

AF  Institute  of  Technology 
Wrlght-Patterson  AFB,  Ohio 
1  Attn:  AFIT  (Library) 

Executive  Director 
AF  Office  of  Scientific  Res. 
Washington  25,  D.C. 

1  Attn:  SREE 

Am  (WLL) 

2  Kirtland  APB,  New  Mexico 
D1 recto 

Air  University  Library 
Maxwell  AFB,  Ala. 

1  Attn:  CR-4582 
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Couunder,  AF  Cambridge  Rea.  Labs. 
ARDC,  L,  0.  Hanacom  Field 
Bedford,  Maaa. 

1  Attn:  CRTOTT-2,  Electronica 


Director 

National  Security  Agency 
Fort  George  G.  Meade,  Md. 
1  Attn:  R42 


Columbia  University 
Dept,  of  Electrical  EngineerinR 
New  York.  N.Y.  10027 
1  Attn:  Prof.  Ralph  J .  Schwarz 


Hqa.,  AF  Byatems  Command 
Andrews  AFB 
Washington  25,  D.C. 

1  Attn:  SCTTAE 


NASA,  Goddard  Space  Flight  Center 
Oreenbelt,  Md. 

1  Attn:  Code  611,  Dr.  G.  H.  Ludwig 


Cornell  U 

Cognitive  Systems  Research  Program 
Ithaca,  N.Y. 


Chief,  Data  Systems  Divisions  1  Attn:  F.  Rosenblatt.  Hollister  Hall 


Asst.  Secy,  of  Defense  (R  and  D) 

R  and  D.  Board,  Dept,  of  Defense 
Washington  25,  D.C. 

1  Attn:  Tech.  Library 

Office  of  Director  of  Defense 
Dept,  of  Defense 
Washington  25,  D.C. 

1  Attn:  Research  and  Engineering 

Institute  for  Defense  Analyses 
400  Army-Navy  Dr.  t 
Arlington,  Va  . ,  222< 

1  Attn:  W.  E.  Bradley 

Defense  Communications 
Dept,  of  Defense 
Washington  25,  D.C. 

1  Attn:  Code  121A,  Tech.  Lib: 

Advisory  Group  on  Electron  De  ^ces 
346  Broadway,  8th  Floor  East 
New  York  13,  N.Y. 

2  Attn:  H.  Sullivan 

Advisory  Group  on  Reliability  of 
Electronic  Equipment 
Office  Asst.  Secy,  of  Defense 
The  Pentagon 

1  Washington  25,  D.C. 

Commanding  Officer 
Diamond  Ordnance  Fuze  Labs. 
Washington  25,  D.C. 

2  Attn:  ORDTL  930,  Dr.  R.  T.  \i 

Diamond  Ordnance  Fuze  Lab. 

U.S.  Ordnance  Corps 
Wasnington  25,  D.C. 

1  Attn;  ORDTL-430-638 

Mr.  R.  H.  Comyn 

Director  ,! 

Weapons  Systems  Evaluation  Group 
1  Washington,  D.C.  ^0305 

U.S.  Dept,  of  Commerce 
National  Bureau  of  Standards 
Boulder  Labs 

Central  Radio  Propagation  Lab. 

1  Boulder,  Colorado 

2  Attn:  Miss  J.V.  Lincoln,  Chief 

RWS8 

NSF,  Engineering  Section 
1  Washington,  D.C. 

Information  Retrieval  Section 
Federal  Aviation  Agency 
Washington,  D.C. 

1  Attn:  MS-112,  Library  Branch 

DDC 

Cameron  Station 
Alexandria  4 ,  Va . 

2  Attn:  TISIA 

U.S,  Coast  Guard 
1300  E.  Street,  N.W. 

Washington  25,  D.C. 

1  Attn:  EEE  Station  5-5 

Chief,  Input  Section 
Clearinghouse  for  Federal  Scientific 
and  Technical  Information,  CFSTI 
SILLS  Building 
5285  Port  Royal  Road 
1  Springfield,  Virginia  22151 


NASA 

Office  of  Adv.  Res.  and  Tech. 
Federal  Office  Bldg.  10-D 
600  Independence  Ave. 

Washington,  D.C. 

1  Attn:  Mr.  Paul  Johnson 

Chief,  U.S.  Army  Security  Agency 
Arlington  Hall  Station 

2  Arlington  12,  Virginia 

Director 

Advanced  Res^rch  Projects  Agency 
1  Washington D.C. 


Asyjjmihsl 

L  if'tx.-  iF  i,i 


^  or  l-LUm^  rkqumlilB 
Ygpp^inr  ■.  ^'ensdi 
Ailk-  HF  4  V.  7Ii-^4W 

llforula  Institute  of  Technology 
PaWdena ,  Calif. 

1  AttX  Prof.  R.  W.  Gould 
1  \  Braverman,  EE  Dept. 

Callfon^  Institute  of  Technology 
4800  Oak\rove  Drive 
Pasadena  ^Callf. 

I  Attn:  Llbi\fy,  Jet  Propulsion  Lab. 


U.  of  Calif  on 
Berkeley  4,  CaN) 

1  Attn;  Prof.  R.lK  luiinders 
Dr.  R.K.  ti>**Tllng, 
Radiation  iJV,  Info.  Dlv. 
Bldg.  30 


Thayer  School  of  Engr. 

Dartmouth  College 
Hanover,  New  Hampshire 
1  Attn:  John  W.  Strohbchn 
Asst.  Professor 

Drexel  Institute  of  Technology 
Philadelphia  4,  Pa. 

1  Attn:  F.  B.  Haynes,  EE  Dept. 

U  of  Florida 

Engineering  Bldg.,  Rm.  336 
Gainesville,  Fin, 

1  Attn;  M.J.  Wiggins,  EE  Dept. 

Georgia  Institute  of  Technology 
Atlanta  13,  Gb  . 

1  Attn:  Mrs.  J.H.  Crosland,  Librarian 
1  F.  Dixon,  Engr.  Experiment 

Station 

Harvard  U 
Pierce  Hall 
Cambridge  38,  Mass. 

1  Attn:  Dean  H.  Brooks,  Dlv.  of  Engr. 

and  Applied  Physics,  Rm.  217 

2  E.  Farkas,  Librarian,  Rm.  303A , 
Tech.  Reports  Collection 

U  of  Hawaii 
Honolulu  14,  Hawaii 
1  Attn:  Asst .  Prol  .  K. 


Nil  )itn  ,  KE  Dept  . 


U  of  minolK 
Urbnno ,  lllinolH 

1  Attn:  P.D.  Coleman,  KR  Res,  Lab, 

1  W.  Perkins,  EE  Res.  Liih. 

1  A.  Albert,  Tcth.Ed..  EE  lie.s. 

1  Library  Serials  . 

1  Prof.  If.  Alpei‘1  ,  Coordinated 

Scl  .  Lab . 


.State  Ualverstlv  ol  Iowa 
Dept,  ol  Electri<-;tl  Engiiiecrl 
lowfi  City,  lo\ia 

EE  Dept.  1  Attn;  Prof.  Dotuild  L.  Epley 


U  of  California 
Los  Angeles  24,  Calif. 

1  Attn:  C.  T.  Leondes,  Prof.  oT 

Engineering,  Engineering  Bept. 
1  R.  S.  Elliot,  Electromagnetles 

Dlv.,  College  of  Engineering 

U  of  California,  San  Diego 
School  of  Science  and  Engineering 
La  Jolla,  Calif. 

1  Attn:  Physics  Dept. 

Carnegie  Institute  of  Tech. 

Schenley  Par.; 

Pittsburg  13,  Pa. 

1  Attn:  Dr.  E.M.  Williams,  EE  Dept. 

Case  Institute  of  Technology 
Engineering  Design  Center 
Cleveland  6,  Ohio 

1  Attn;  Dr,  J.  B,  Neswlck,  Director 

Columbia  Radiation  Lab. 

Columbia  University 
538  W.  120th  St. 

New  York,  N.Y.  10027 

♦No  AF  or  Classified  Reports 


♦Instltuto  tic  Pcsqiiis.i.s  dji  MiU'inha 
MlniBterlo  da  Marinha 
Rio  de>  Janeiro 
Estado  dn  Guanabara,  Brazil 
1  Attn:  lloberto  U.  d;i  Costa 

Johns  Hopkln.s  U 
Charles  and  34th  St. 

Baltimore  IK,  Md, 

1  Attn:  Librarian,  ('iirlvlc  Barton  Lab. 

Johns  Hopklii.s  V 
8621  Georgia  Avo. 

Silver  Sprltiga,  Md , 

1  Attn:  N.  H.  rhoks> 

1  Mr.  A, Vi.  Si'gy,  Applied 

Physlc‘3  Lab, 

1  Attn;  Mr.  E.  E.  Green 

Llnfleld  Hesoarch  Institulo 
McMinnvt  I  le  ,  Ore . 

1  Attn;  G.  N.  Hlekok,  Dlrec  lor 


Murtiuottu  Uni\t*istt\ 

College  of  ling  1  iieer  1 
1515  W.  Wiseotisln  A\’c-. 
Milwaukee  3,  Wls, 

1  At  tn;  A  .C.  Mod  ler,  Lb  [fept  . 
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MIT 

Cunbrldge  39 »  Maas. 

1  Attn:  Raa.  Lab.  of  Elec.,  Doc.  Rn 
1  Mlaa  A.  Slla.  Llbn.  Rm.  4-244, 

LIR 

1  Mr.  J.E.  Ward,  Elac.Sya.  Lab. 

H  1  T 

Lincoln  Laboratory 
P.O.  Box  73 

1  Attn:  Lexington  73,  Masa. 

1  Navy  Repraaentatlve 

1  Dr.  W.I.  Wells 

1  Kenneth  L.  Jordon,  Jr. 


«HJ  of  Saskatchewan 
College  of  Engineering 
Saakatooi,  Canada 
1  Attn:  Prof.  R.  £.  Ludwig 

Syracuao  U 
Syracuse  10.  N.Y. 

1  Attn:  EE  Dept. 

MI  of  Texas 
Serials  Acquisitions 
Austin,  Texas 
1  Attn:  Mr.  John  Womack 

Serials  Librarian 


Bell  Telephone  Labs.,  inc. 

Technical  Information  Library 
Whlppany,  N.J. 

1  Attn:  Tech.  Repts.  Llbrn., 

Whlppany  Lab. 

The  Boeing  Company 
Mall  Stop  MS-1331-ORO,  1-8000 
Seattle  24,  Washington 
1  Attn;  Dr,  Ervin  J.  Nalos 

•Central  Electronics  Engineering 
Research  Institute 
Pllanl,  Rajasthan,  India 
1  Attn:  Om  P.  Gandhi  -  Via;  ONR/London 


U.  of  Michigan 
Ann  Arbor,  Mich. 

1  Attn:  Dir.,  Cooley  Elec  Labs,, 
N,  Campus 


•Uppsala  U 

Institute  of  Physics 
Uppsala,  Sweden 
1  Attn:  Dr.  P.  A.  Tove 


U  of  Minnesota 
Institute  of  Technology 
Minneapolis  14,  Minn. 

1  Attn:  Prof.  A.  Van  der 
EE  Dept. 


1  Dr.  J.E.  Rowe,  Elec.  Phys.  Lab. 

1  Comm.  Scl.  Lab,,  180  Prleae  Bldg. 

U  of  Michigan 

Institute  ofs^clence  and  Technology  1 
P.O.  Box  618 
Ann  Arbor,  Mich 
1  Attn:  Tech. 

I  W. 


U  of  Toledo 

Dept,  of  Electr.  Engr. 

Toledo  6,  Ohio 

Attn:  James  B.  Farlson 


U  of  Nevada 

College  of  Engineering 
Reno,  Nevada 
1  Attn:  Dr.  R.A.  Manhart,  EE  Dept 

Northeastern  U 
The  Dodge  Library 
Boston  15.  Mess. 

1  Attn:  Joyce  E.  Lunde,  LIhririM 


Northwestern  U 
2422  Oakton  St . 

Evanston,  III. 

1  Attn:  W.S.  Toth  Aerial 
Measurements  Lab. 


U  of  Notre  Dame 
South  Bend,  Ind. 

1  Attn;  E.  Henry,  EE 


Ldhiary 


Ohio  State  U 
2024  Niel  Ave. 

Columbus  10,  Ohio 

1  Attn:  Prof.  E.M.  Boone,  EE  Dept. 


Oregon  State  U 
Corvallis,  Ore. 

1  Attn:  H.J.  Oorthuys,  EE  Dept. 


ILV  ■  I.UH 

Ibl 

g  itta:  Mr  d.t 


Airborne  Instruments  Lab. 


U  of  Notre  Dame 
Dept,  of  Electrical 
Notre  Dame,  Indiana 
1  Attn:  W.  B.  Berry 


Admiral  Corp. 

3800  Cortland  St. 
Chicago  47 ,  III . 

1  Attn:  E.N.  Robertson, 


Polytechnic  Institute 
333  Jay  St . 

Brooklyn,  N.Y. 

1  Attn:  L.  Shaw,  EE  Dept. 


Comae  Road 

Deer  Park.  Long  Island,  N.Y. 

1  Attn:  J.  Dyer,  Vice-Pres.  and 
Tech .  Dir. 


Polytechnic  Institute  of  Brooklyn 
Graduate  Center,  Route  110 
Farmingdale,  N.Y. 

1  Attn:  Librarian 


Autonetlcs 

Dlv.  ol  North  American  Aviation,  Inc. 
9150  E,  Imperial  Highway 
Dfwney,  Calif. 

1  Attn:  Tech.  Library  3040-3 


1 


1 


1 


1 


1 


1 


1 


1 


1 

1 


1 


1 


I 


Convalr  -  San  Diego 

Dlv.  of  General  Dynamics  Corp. 

San  Diego  12,  California 
Attn:  Engineering  Library 

Cook  Research  Labs. 

6401  W.  Oakton  St. 

Attn:  Morton  Grove,  Ill. 

Cornell  Aeronautical  Labs.,  Inc. 
4455  Genessee 
Buffalo  21,  N.Y. 

Attn:  Library 

Eltel-McCullough,  Inc. 

301  Industrial  Way 
San  Carlos,  Calif. 

Attn:  Research  Librarian 

Ewan  Knight  Corp. 

East  Natick.  Mass. 

Attn:  Library 

Fairchild  Semiconductor  Corp. 

313  Fairchild  Dr. 

P.O.  Box  880 
Mt.  View.  Calif. 

Attn:  Dr.  V.  H.  Orlnlch 

General  Electric  Co. 

Defense  Electronics  Dlv.,  IMED 
Cornell  University,  Ithaca,  N.Y. 
Attn:  Library 

Via:  Commander,  ASD  W-P  AFB,  Ohio 
ASRNGW  D.E.  Lewis 

General  Electric  TWT  Products  Sec. 
601  California  Ave. 

Palo  Alto,  Calif. 

Attn;  Tech.  Library,  C.  0.  Lob 

General  Electric  Co.  Res.  Lab. 

P.O.  Box  1088 
Schnectady,  N.Y. 

Attn;  Or.  P.M.  Lewis 

R.  L,  Shuey,  Mgr.  Info. 
Studies  Sec . 

General  Electric  Co. 

Electronics  Park 
Bldg.  3,  Rm.  143-1 
Syi’acuse,  N.Y. 

Attn:  Doc.  Library,  Y.  Burke 

Gllfillan  Brothers 
1815  Venice  Blvd. 

Los  Angeles,  Calif. 

Attn:  Engr.  Library 

The  Halllcraf tors  Co. 

5th  and  Kostner  Ave. 

All  II :  Chicago  24  ,  III  . 


Purdue  U 
Lafayette,  Ind. 

1  Attn:  Library,  EE  Dept. 

Rensselaer  Polytechnic  Institute 
School  of  Engineering 
Troy,  N.Y. 

1  Attn:  Library,  Serials  Dept. 

Kenneth  E.  Uortenson 


Dell  Telephone  Labs. 
Murray  Hill  I.nb. 


Murray 

Hill , 

N.J. 

1 

Attn: 

Dr. 

.1 . 

, R.  Pierre 

1 

I)r. 

S, 

.  Darlington 

1 

Mr. 

A 

.  .1.  Grossman 

•No  AF  or  Classified  Reports. 


Hewlett-Packard  Co. 

1501  Page  Mill  Hoad 
1  Attn;  Palo  Alio,  Calif. 

Hughes  Alrcrait 
Malibu  Beach,  Calif. 

1  Attn:  Mr .  lams 
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Hughes  Aircraft  Co. 

Florence  at  Teale  St. 

Culver  City,  Calif. 

1  Attn;  Tech.  Doc.  Cen . ,  Bldg.  6 
Rm  C2046 


*Dlr.,  National  Physical  Lab. 
Hilslde  Road 
New  Delhi  12,  India 
1  Attn;  S.C.  Sharma  -  Via: 

ONR /London 


Sylvania  Electronics  System 
100  First  St. 

Waltham  54.  Mass. 

1  Attn;  Librarian,  Waltham  Labs 
I  Mr.  E.E.  Hollis 


Hughes  Aircraft  Co. 

P.O.  Box  278 
Newport  Beach,  Calif. 

1  Attn;  Library,  Semiconductor  Div. 


♦Northern  Electric  Co.,  Ltd. 
Research  and  Development  Labs. 
Dept .  8300 
Ottawa,  Ontario 
1  Attn:  H.  L.  Blacker 


IBM,  Box  390  Boardman  Road 
Poughkeepsie,  N.Y. 

1  Attn:  J.C.  Logue,  Data  Systems  Div, 

IBM  Poughkeepsie,  N.Y. 

1  Attn:  Product  De’*.  bau., 

E .  M .  Da  h'i  8 

IBM  ASD  and  Research  Library 
P.O.  Box  66 

Los  Gatos,  Calif  95031 
1  Attn:  Miss  K.  Griffin,  Bldg.  025 

ITT  Federal  Labs. 

500  Washlngt>m  Ave 
Nutley  10,  N.-T. 

1  Attn:  Mr,  K.  Mount,  Ll'IirBrlan 

Laboratory  for  Electronici 
1075  Commonwealth  Ave. 

Boston  15,  Mass. 

1  Attn:  Library 

LEL,  Inc. 

75  Akron  St , 

Coplague,  Long  Island,  N.Y 
1  Attn:  Mr.  R.  S.  Mautner 


Lenkurt  Electric  Co. 

San  Carlos,  Calif. 

1  Attn:  M.L.  Waller.  Librarian 

Llbrascope 
Div.  of  General  Precision,  I 
806  Western  Ave. 

Glendale  1,  Calif. 

1  Attn:  Engr.  Library 


Lockheed  Missiles  and  Spa 
P.O  Box  304,  Bldg.  524 
Sunnyvale,  Calif. 

Attn:  Dr.  W.M.  Harris, 

0.  w.  Price,  De^ 

Melpar,  Inc.  / 

3000  Arlington  Blvd./ 
Falls  Church,  Va.  / 
Attn:  Librarian  / 


Microwave  Associates,  Inc. 
Northwest  Industrial  Park 
Burlington,  Mass. 

1  Attn;  K.  Mortenson 
1  Librarian 

Microwave  Electronics  Corp. 
4061  Transport  St. 

Palo  Alto,  Calif. 

1  Attn;  S.F.  Kalsel 
M.C.  Long 


Northronlcs 

Palos  Verdes  Research  Park 
6101  Crest  Ro 

Palos  Fnrita^^ltates,  Calif. 

1  Attn:  l^o.  Center 

Pacific -r^'isrkiductors,  Inc. 

14520  Af  Li  tlon  Blvd. 

Lawndal/,  E^|.if. 

1  Attn:  til  PUrth 

/ 

Mihi^o,  Tech.  Rep.  Division 
>,IU  Box  10 

Vashlngton,  Pa. 

I  l/l:n:  F.  R.  Sherman 

ykadlo  Corp.  of  Amerl  :a 
RCA  Labs.,  David  Sartoff  Res.  Cen. 
Princeton,  N.J. 

Attn:  Dr.  J.  Sklansky 

■a  Lslim  ,  -PnimaTaii,  H.4, 
t  Am  M  iMIisfii 


■j-ar.  Librarian 


Raytheon  Manufacturing  C&v 
Microwave  and  Power  Tube  Div. 
Burlington,  Mass. 

1  Attn:  Librarian,  Spencer  Lab. 

Raytheon  Manufacturing  Co. 

Res.  Div.,  28  Seyon  St. 

Waltham,  Mass. 

1  Attn:  Dr.  H.  Statz 
1  Mrs.  M.  Bennett,  Librarian 

1  Research  Div.  Library 

Sandla  Corp. 

Sandla  Base,  Albuquerque,  N.M. 

1  Attn;  Mrs.  B.  R.  Allen,  Librarian 

Scientific  Atlanta.  Inc. 

P.O.  Box  13654 
Atlanta,  Georgia  30524 
1  Attn:  Dr.  John  E.  Pippin 

Director  of  Research 

Sperry  Rand  Corp. 

Sperry  Electron  Tube  Div. 
Gainesville,  Fla. 

1  Attn:  Librarian 


Mlnneapolls-Honeywell  Regulator  Co. 
1177  Blue  Heron  Blvd. 

Riviera  Beach,  Fla, 

1  Attn:  Semiconductor  Products  Library 

The  Mitre  Corp. 

Bedford,  Mass. 

1  Attn:  Library 

Monsanto  Research  Corp. 

Station  B,  Box  8 
Dayton  7,  Ohio 
1  Attn:  Mrs.  D.  Crabtree 

Monsanto  Chemical  Co. 

800  N.  Lindbergh  Blvd. 

St .  Louis  66,  Mo. 

1  Attn:  Mr.  E.  Orban,  Mgr. 

Inorganic  Dev, 


Sperry  Gyroscope  Co. 

Div,  of  Sperry  Rand  Corp, 

Great  Neck.  N.Y. 

1  Attn:  L.  Swern  (MS3T10S) 

Sperry  Gyroscope  Co. 

Engineering  Library 

Mall  Station  F-7 

Great  Neck,  Long  Island,  N.Y. 

1  Attn:  K.  Barney,  Engr.  Dept.  Head 

Sylvania  Electric  Products,  Inc. 
500  Evelyn  Ave. 

1  Mt.  View,  Calif. 

1  Attn:  Mr,  E.  0.  Ammann 


♦No  AF  or  Classified  Reports. 


Technicnl  Resenrch  Group 
Route  No.  110 

1  Melville,  New  York  11749 

Texns  Instruments,  Inc. 

P.O,  Box  6015 
Dallas  22,  Texas 

1  Attn:  M.E.  Chun,  Apparatus  Div. 

Texas  Instruments,  Inc. 

P.O.  Box  3012 
Dallas.  Texas  75222 

1  Attn:  Tech.  Repts.  Service,  MS-65 

2  Semi-Conductor  Components 
Library 

Texas  Instnjments,  Inc. 

Corporate  Research  and  Engineering 
Technical  Reporl.s  Service 
P.O.  Box  5474 
1  Dallas  22,  Texas 

Tektronix,  Inc, 

P.O.  Box  500 
Beaverton,  Ore. 

4  Attn:  Dr,  J.F.  DeLord,  Dir.  of 
Research 

Varlan  Associates 
611  Hansen  Way 
Palo  Alto,  Calif. 

1  Attn:  Tech.  Library 

Weltermann  Electronics 
4549  North  3Hth  St. 

1  Milwaukee  9,  Wisconsin 

Westlnghousc  Electric  Corp. 
Friendship  International  Airport 
Box  746,  Baltimore  3,  Md . 

1  Attn:  G.R,  Kilgore,  Mgr.  Appl .  Res. 
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